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ABSTRACT 

Analysis  of  survival  constraints  on  the  design  of 
solid  targets  for  tokamak  bundle  divertors  is  presented. 
Previous  target  design  efforts  are  reviewed.  Consider¬ 
ations  of  heat  removal,  surface  erosion,  and  fatigue 
life  are  included  in  a  generalized  design  window 
methodology  which  facilitates  target  selection.  Using 
subcooled  water  as  coolant,  eight  possible  target 
materials  are  evaluated  for  use  in  tubular  and  plate 
targets  as  substrates,  coatings,  and  claddings.  Subject 
to  the  severe  environment  of  the  tokamak  plasma,  the 
most  promising  conventional  designs  are  identified. 

A  thermally  bonded,  mechanically  unbonded  laminated 
design  is  proposed  and  evaluated  as  a  target  design 
well  suited  to  the  divertor  target  environment.  Due 
to  fatigue  and  sputtering  erosion  this  configuration 
has  limited  life,  but  appears  to  constitute  an  upper  . 
bound  for  the  capabilities  of  a  solid  target  design. 
Needs  for  experimental  work  are  identified. 


Thesis  Supervisor:  Prof.  Neil  E.  Todreas 


,  -cession  For 

rr*G  gka&i 
j  r  ••••'  t'.b 

i  -m 

i  :  v.icn. 


n 

o 


ACKNOWLEDGEMENT 


The  author  appreciates  the  patience  and  assistance 
of  the  many  who  have  helped  make  this  project  possible. 

In  particular,  the  wisdom  and  suggestions  of  my  advisors. 
Professor  B.B.  Mikic  and  Professor  N.E.  Todreas,  were 
priceless  aids.  The  vigor,  insight,  and  helpfulness  of 
Paul  Gierszewski  knows  no  equal.  Ted  Yang  and  Bruce 
Montgomery  offered  support  and  key  direction  in  the 
formulation  of  this  project. 


4CJk 


TABLE  OF  CONTENTS 


9  Title  Page 
Abstract 
Acknowledgement 
Table  of  Contents 
Nomenclature 

Introduction 

A.  Functions  of  Divertor 

B.  Focus  of  Work 

C.  Summary  of  Other  Divertor 

Work  to  Date 

D.  Base  Case  Parameters 

Thermo-Hydraul ic s 

A.  Coolant  Survey 

B.  Determination  of  Applicable 

CHF  Correlations 

C.  Pressure  Drop  Limitations 

D.  Heat  Transfer  Coefficient 

E.  Swirl  and  Mixed  Flow  Schemes 

F.  Effect  of  Asymmetric  Heating 

Surface  Interactions 

A.  Sputtering  and  Plasma 

Contamination 

B.  Internal  Erosion 

Thermo-Mechanics 

A.  Pressure  Stress 

B.  Thermal  Stress 

C.  Fatigue 

D.  Thermal  Shock 

E.  Temperature  Limits 

F.  Graphical  Materials  Data 

Evaluation  of  Design  Options 

A.  Possible  Geometries 

B.  Materials 

1.  Functional  Requirements 

2.  Limits  of  Data  Base 

C.  Particle  Load  vs.  Lifetime 

Tradeoff 

D.  Single  Material  Analysis 

1.  Tube  Design 

2.  Plate  Design 

E.  Composite  Material  Analysis 

1.  Co "  *-  -  d  Tube  C-; 

2.  Armored  Plate  uesign 


77 


3.  Mechanically  Unbonded  Layered 
Plate  Design 


Sensitivity  Analysis 

92 

"  A. 

Pulse  Length  and  Fatigue 

Constraint 

92 

B. 

Sputtering  Characteristics 

93 

C. 

CHF  Limit 

94 

D. 

Pumping  Power  Limit 

95 

E. 

Inlet  Temperature 

96 

F. 

Pressure 

96 

G. 

Tube  Length 

97 

H. 

Thickness 

97 

I. 

Heat  Flux 

98 

Conclusions 

109 

Recommendations 

111 

References 

Appendix  A:  Materials  Properties 
Appendix  B:  Methodology  Summary 


113 

118 

121 


NOMENCLATURE 


at  * 

a 

Thermal  diffusivity 

s  . 

mt 

Allowable  stress 

Tube  inner  radius 

intensity 

b 

Tube  outer  radius 

S.F. 

CHF  safety  factor 

c_ 

Specific  heat 

t 

Wall  thickness 

P 

C 

Target  capacity 

T 

Temperature 

factor  (fraction  of 

Ti 

T. 

Ion  Temperature 

year  exposed  to 

Fluid  inlet 

plasma) 

in 

Temperature 

D 

Channel  hydraulic 
diameter 

4Tl-2 

Temperature  rise 
in  coolant  from 

E 

Modulus  of  elasticity 

inlet  to  outlet 

f 

Friction  factor 

V 

Coolant  velocity 

g„ 

Gravitation  constant 

V* 

Pumping  power  to 

Mass  flux 

p 

heat  transfer 

G 

ratio 

h 

Heat  transfer 
coefficient 

y 

Tape  twist  ratio 
(inside  dia./ 

hfg 

Enthalpy  of 

180°  twist) 

evaporation 

oc 

Thermal  expansion 

x 

Max  allowable  impur- 

coefficient 

max 

ity  in  plasma  (%) 

c 

Strain 

Particle  flux  for  i'th 
isotope 

e 

Angle  from  normal 

k 

Thermal  conductivity 

M 

Viscosity 

kn 

Boltzmann's  constant 

D 

Poisson's  ratio 

B 

L 

Channel  length 

f 

Density 

Mass 

<r 

Stress 

m 

• 

m 

Mass  flow  rate 

r 

Time 

M 

Atomic  weight 

U) 

Period  of  heat 
load  oscilla- 

N 

Particle  current 

tion 

na 

Avogadro ' s  Number 

Subscripts : 

P 

Pressure 

b 

Bulk  fluid 

P 

Power 

c 

Critical 

q" 

Heat  flux 

s 

Surface 

r 

Channel  inner  radius 

w 

Fluid  at  wall 

S. 

Sputtering  coefficient 

temperature 

1 

(atom/ion) 

I.  INTRODUCTION 


A.  Functions  of  the  Divertor 

In  order  to  advance  toward  eventual  commercial  power 
output  from  tokamak  fusion  machines,  methods  must  be  devised 
to  permit  uninterrupted  operation  of  devices  over  long 
periods  of  time.  Only  in  this  way  can  reliable  and  steady 
electrical  energy  be  supplied  to  the  distribution  grid. 
Current  experimental  machines  are  capable  of  fusion  pulses 
of  durations  less  than  one  second.  Pulse  duration  is 
presently  limited  by  the  flux  swing  in  the  magnet  power 
supply.  If  this  obstacle  can  be  overcome,  two  secondary 
limits  on  extending  the  pulse  duration  result  from 
(a)  the  generation  by  the  fusion  process  of  helium, 
which  quenches  the  plasma  burn  by  displacing  hydrogen 
fuel,  and  from  (b)  the  sputtering  of  impurity  atoms  from 
the  chamber  first  wall,  which  quench  the  plasma  by 
greatly  enhanced  line  radiation.  An  additional  related 
constraint  is  the  damage  sustained  by  the  first  wall 
itself  due  to  sputtering,  necessitating  frequent  and 
costly  refurbishment. 

The  development  of  the  divertor  provides  a  potential 
remedy  for  these  problems.  In  general,  a  divertor  must 
serve  as  both  a  "helium  exhaust"  system  and  an  "impurity" 
filter"  system  for  the  tokamak  plasma.  An  outer  sheath 

Ou  the  toroidal  plasma  column  is  magnetically  diverted 


Scrape-off  Layer 


Mo)or  radius  (m) 


Figure  I— 1 :  Horizontal  section  of  bundle  divertor  showing 
two  additional  exhaust  coils  to  expand  the 
diverted  flux  bundle,  *■) 


to  a  separata  chamber,  where  these  vital  particle 

handling  functions  may  be  undertaken  by  impingement  on  an 

5 

appropriate  target.  (See  Pigs.  I-l,I-2.)  Helium  is 
continuously  collected  as  it  diffuses  outward  from  the 
reacting  plasma  core;  impurities  are  filtered  before  they 
can  penetrate  to  the  energetic  core,  and  the  rate  of 
damage  to  the  chamber  first  wall  is  significantly  reduced. 
The  reason  for  a  divertor's  existence,  then,  is  this 
particle  handling  function.  Numerous  methods  have  been 
proposed  to  carry  out  this  task,  including  chemical  and 
mechanical  gettering,  diffusion  of  light  atoms  through 
metals,  and  neutral  gas  pumping.  Each  method  implies  use 
of  some  target  machanism  which  intercepts  the  diverted 
plasma . 

The  environment  in  which  the  divertor  target  must 

perform  will  be  extremely  severe.  Though  estimates  vary, 

2 

heat  loads  on  the  order  of  1  kw/cm  ,  particle  loads  of 
22  2 

10  particles/m  -sec,  and  particle  energies  from  1  to 
10  keV  are  anticipated  to  impinge  on  the  target.  In  order 
to  permit  economical  long  term  operation  of  the  fusion 
device,  the  system  must  be  designed  to  avoid  excessive 
down  times  for  target  refurbishment.  The  McDonnell- 
Douglas  study  of  tokamak  maintainability  suggests  a  min¬ 
imum  period  between  maintenance  shutdowns  of  one  year  as 
essential  to  fusionTs  economic  viability. (Ref.  2)  Thus, 
a  prerequisite  to  rilling  i — :  particle  handling 


function  is  the  ability  to  survive  the  extreme  environment 

for  at  least  a  one  year  period.  To  facilitate  development 

v 

of  a  workable  divertor  target,  designs  which  appear  to 
have  adequate  longevity  will  be  sought.  The  particle 
handling  function  can  then  be  added  to  a  sufficiently 
durable  target.  This  decoupling  of  energy  absorption 
from  particle  handling  may  prove  attractive  in  future 
designs. 

B.  Focus  of  Work 

While  possibilities  for  divertor  target  design  in¬ 
clude  novel  concepts  such  as  liquid  targets,  gas  targets, 
and  curtains  of  moving  pellets,  (see  Ref.  3)  the  conven¬ 
tional  design  encompasses  a  solid  target  material  which 
acts  at  a  minimum  to  absorb  energy,  and  possibly  also  to 
trap  particles.  This  paper  will  concentrate  on  various 
arrangements  of  solid  target  materials  under  active  cool¬ 
ing  in  an  effort  to  develop  a  more  comprehensive  knowledge 
of  design  tradeoffs  among  conventional  solid  target 
options.  A  design  window  approach  will  help  define  the 
possible  operating  configurations.  The  major  constraints 
on  survivability  are  heat  removal,  surface  erosion  rates, 
and  fatigue  life.  Known  data  and  materials  behavior  for 
eight  selected  materials  will  be  incorporated  in  developing 
a  design  methodology  for  optimizing  target  life  subject  to 
heat  load  and  pumping  power  constraints,  and  for  optimiz¬ 
ing  the  allowable  heat  load  subject  to  lifetime  constraints. 


The  eight  candidate  elements  and  alloys  are:  Beryllium, 
Graphite,  Aluminum,  Titanium  Alloy  (Ti-6A1-4V) ,  Vanadium 
Alloy j (V-Z5Cr-.th)  ,  Copper  Alloy  (ZAC-2),  Niobium  Alloy  (D- 
43) ,  and  Molybdenum  Alloy  (TZM) . 

In  considering  heat  removal,  several  thermo- 
hydraulic  concerns  must  be  addressed.  The  operating 
limits  imposed  by  critical  heat  flux,  fluid  pressure  drop, 
and  heat  transfer  coefficient  must  be  defined.  The  poten¬ 
tial  benefits  of  swirl  and  mixed  flow  schemes  and  the 
impact  of  axisymmetric  heating  will  be  assessed.  A  survey 
of  possible  coolants  indicates  that  subcooled  v/ater  offers 
advantages  over  other  coolants  in  handling  high  heat 
fluxes.  Water  coolant  boasts  the  broadest  range  of 
correlated  behavior,  and  will  be  used  exclusively  for  this 
study. 

Sputtering  rates  based  on  controlled,  single 
incident-species  experiments  will  be  extrapolated  to  the 
projected  environment  of  the  next  generation  tokamak 
experiments  to  indicate  the  severity  of  this  major  con¬ 
straint  on  the  target  design  process.  The  impact  on  the 
plasma  and  on  particle  handling  requirements  will  also 
be  examined. 

Thermo-mechanics  of  solid  targets  will  be  treated  in 
accordance  with  guidelines  of  the  ASME  Pressure  Vessel 
Code  for  Nuclear  Components.  Tube  and  flat  plate 

rec-Te tries  will  be  ev? lusted  for  cases  of  a  sincle 
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material  ,  a-.d  for  cases  of  a  coating  or  cladding  over  a 
substrate. 

the  divertor  target  design  procedures  developed  here 
should  be  applicable,  with  minor  modification,  to  the 
development  of  actively  cooled  limiters  and  divertor  duct 
walls. 

C.  Summary  of  Other  Divertor  Work  to  Date 

Several  detailed  solid  divertor  target  designs  have 
been  produced  to  date.  These  range  in  intended  use  from 
developmental  programs  with  low  capacity  factors  to 
commercial  power  reactors  with  stringent  availability 
requirements.  Each  is  currently  subject  to  some  doubt  as 
to  survivability  at  the  commercial  scale. 

The  Poloidal  Divertor  Experiment  is  now  operating  at 
Princeton  to  assess  the  feasibility  of  employing  divertors 
as  impurity  and  exhaust  control  devices.  Using  titanium 
plates  as  particle  collectors,  the  experiment  has 
successfully  performed  the  divertor  functions  for  very 
low  energy  plasma  edge  particles  (20  to  30  eV) . 

The  category  of  developmental  systems  includes  the 
proposals  of  Grumman  (SLPX,  Ref.  4),  Oak  Ridge  (ETP,  Ref. 
5),  and  the  U.S.  INTOR  report  (Ref.  6).  The  Grumman  design 
for  the  conceptual  Superconducting  Long  Pulse  Experiment 
at  Princeton  envisions  a  solid  poloidal  divertor  which 
incorporates  swirl  water  flow  in  a  ladder-like  tube  array 
to  cool  a  total  heat  load  of  40  MW.  Pulse  durations  are 


30  seconds.  Average  best  flux  on  the  tubes  is  1.5  kw/cm2. 

Gettering  of  reflected  particles  and  regeneration  heating 

is  th4  projected  particle  handling  method.  The  heat 

removal  system  employs  a  design  previously  demonstrated 

2 

to  handle  a  heat  load  of  3.2  kw/cm  (Ref.  7),  but  no 

consideration  of  surface  sputtering  or  fatigue  failure  is 

made.  The  Oak  Ridge  design  handles  a  total  heat  load  of 

2 

50  MW,  an  average  heat  load  of  1.1  kw/cm  ,  and  a  pulse 

4 

length  of  100  seconds.  It  will  undergo  5  x  10  cycles  per 

year.  The  INTOR  design  proposes  tungsten  protective  tile 

brazed  over  water  cooled  copper  plates  to  accomodate  an 

2 

average  heat  load  of  0.5  kw/cm  .  The  one  centimeter  thick 

tiles  are  projected  to  last  two  years  in  the  sputtering 

environment,  but  the  effect  of  sputtered  tungsten  on  plasma 

energy  is  a  major  unanswered  question. 

In  a  more  generalized  category  of  divertor  target 

design  is  the  proposal  of  Westinghouse  (Ref.  8),  a  bundle 

divertor  target  consisting  of  sloped  tube  arrays  which 

2 

handle  heat  fluxes  up  to  3.2  kw/cm  .  Special  considera-' 

tions  to  protect  tube  bends  exposed  to  cyclic  thermal 

stresses  and  sputtering  have  not  been  provided  for. 

On  the  commercial  scale,  UWMAK  III  advocates  a 

poloidal  divertor  with  energy  absorbing  TZM  collector 

plates.  Sodium  coolant  is  used  to  handle  an  average 

2 

heat  load  of  0.589  kw/cm  .  Burn  length  is  30  minutes. 

Energetic  incident  particles  impinge  on  a  thin  sacrificial 


8L_k. 
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T2M  sheet  which  must  be  replaced  in-situ  every  17  hours. 

Particle  collection  is  then  achieved  by  a  series  of 
f 

cryopumps . 

Facilities  are  in  existence  which  handle  heat  fluxes 
2 

from  2  to  12  kw/cm  for  specific  non-tokamak  configurations, 

notably  those  of  NASA- Ames  (Ref.  9) ,  McDonnell-Douglas 

(Ref.  10),  and  RCA  (Ref.  11).  Fatigue  and  sputtering 

constraints  are  not  major  concerns  for  these  devices. 

Investigations  into  first  wall  cooling  capabilities 

have  been  made,  resulting  in  specific  design  points  for 

high  heat  flux  situations.  Notable  is  the  report  of 

Hoffman,  et  al.  (Ref.  12),  which  predicts  safe  heat  flux 

2 

capabilities  of  3.3  kw/cm  for  thin  walled,  1  cm  inner 

diameter,  1  meter  long  tubes  heated  on  one  side.  Their 

anticipated  pumping  power  to  heat  transfer  ratio  is  3%. 

No  provision  is  made  for  sputtering  or  fatigue. 

Input  from  the  Divertor  Workshop  at  MIT,  April  1980, 

and  the  Plasma  Materials  Workshop  at  Sandia-Albuquerque, 

June  1980,  led  to  selection  of  a  design  heat  load  of 
2 

1  kw/cm  as  reasonable  for  the  divertor  target  plate,  based 

on  existing  materials  constraints.  Off  normal  conditions 

2 

resulting  in  temporary  heat  loads  up  to  2  kw/cm  were 

deemed  possible.  Experimentation  is  today  under  way  at 

Sandia-Albuquerque  to  achieve  acceptable  operation  under 
2 

a  2  kw/cm  heat  load. 


/ 


The  recent  STARFIRE  proposal  of  Argonne  National 

Laboratory  (Ref.  13)  is  of  interest  to  divertor  designers, 

even  though  the  proposed  commercial-scale  reactor  employs 

a  toroidally  mounted  limiter  to  fulfill  the  divertor 

functions.  Survival  of  the  particle  collecting  plates  is 

greatly  enhanced  by  extracting  a  greater  fraction  of  the 

plasma  energy  as  radiation  and  by  requiring  a  relatively 

low  particle  load  on  the  surface.  The  anticipated  peak 

2 

heat  loads  of  0.4  kw/cm  make  survivability  quite  possible. 
Continuous  operation  also  enhances  survival  chances  by 
greatly  reducing  fatigue  constraints.  Beryllium  coating 
material  ostensibly  overcomes  the  sputtering  erosion 
problem  by  redepositing  beryllium  on  the  limiter.  Each 
of  these  features  is  of  great  value  in  prolonging  divertor 
life. 

D.  Base  Case  Parameters 

As  noted  earlier,  the  desired  divertor  target  in  this 

investigation  must  survive  for  one  year  under  operating 

2 

heat  loads  of  1  kw/cm  .  In  order  to  permit  consistent 

analysis,  operating  conditions  derived  from  ETF  arid  INTOR 

studies  are  established  as  follows. 

-Pulse  duration  is  90  seconds. 

-Rejuvenation  period  is  15  seconds. 

-Three  targets  are  assumed,  such  that  each  spends  1/3 
of  the  time  "on  line"  and  2/3  of  the  time  in  a 
"particle  regeneration"  status. 

-Devices  operate  with  no  prolonged  shutdowns  for  a 


one  year  period. 

5 

-Each  target  undergoes  10  plasma  pulse  cycles  per 
year. 

-Ion  temperature  at  the  target  is  1.3  keV. 

-Particle  current  is  1.5  x  1022  /sec. 

-Available  surface  area  for  a  bundle  divertor  is  9  m  . 

-Particle  flux  is  1.67  x  1022  /sec-m2. 

To  establish  the  design  window  procedure,  a  set  of 
base  case  parameters  for  the  target  is  established.  The 
selection  of  these  parameters  is  guided  by  analysis  of  the 
impacts  of  ion  sputtering,  cyclic  fatigue,  subcooled  boil¬ 
ing  pressure  drop,  and  critical  heat  flux.  They  are 
selected  as  representative  of  the  practical  ranges  of  each 
variable.  The  10  cm  length  is  selected  as  reasonable  for 
near-term  single  component  testing.  Variation  of  each 
parameter  is  examined  in  Chapter  VI.  The  base  case 
parameters  are : 

Heat  Flux  q"  =  1  kw/cm2 

Target  Life  (Goal)  T  -  1  year 

Coolant  Channel  Length  L  -  10  cm 

Channel  Diameter  D  *  1  cm 

Coolant  Pressure  p  **  500  psia, 

3.45  MPa 

Coolant  Inlet  Temoerature  T.  =  30  #C,  86  °F 

ir* 

t+r 

r 


Tube  Outer  Radius  to 
Inner  Radius  Ratio 


1.2 
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II.  THERMO-HYDRAULICS 

A.  Cdolant  Survey 

Successful  heat  removal  for  the  base  case  target 
design  requires  an  actively  cooled  system.  A  survey  of 
a  variety  of  basic  options  for  active  cooling  was  made. 
Liquid  metals,  helium  gas,  subcooled  water,  and  mist-steam 
flows  were  compared  for  viability.  High  pressure  drops 
occur  for  liquid  metals  flowing  transverse  to  the  diverted 
magnetic  field  lines  (Refs.  14,15).  Helium  offers  safety 
and  materials  compatibility  advantages,  but  was  assessed  as 
inferior  to  subcooled  water  in  the  extensive  analysis  of 
the  STARFIRE  study  group  (Ref.  13) .  Mist-steam  flow  can 
operate  at  lower  coolant  pressures  than  is  practical  for 
subcooled  water,  but  materials  are  not  known  which  could 
operate  at  the  high  resulting  wall  temperatures.  Mist- 
steam  flows  would  also  demand  a  relatively  high  pumping 
power.  Subcooled  water  emerges  as  the  most  attractive 
coolant  for  the  high  heat  flux  divertor  target,  assuming 
adequate  provision  for  safeguarding  against  water-alkali 
metal  reactions  is  possible. 

B.  Determination  of  Applicable  CHF  Correlations 

Under  the  very  high  heat  loads  experienced  by  the 
divertor  target,  channels  cooled  by  subcooled  liquid  are 
subject  to  catastrophic  failure  resulting  from  coolant 
burnout.  Accordingly,  the  first  limitation  on  fluid 
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operating  parameters  is  che  proximity  of  the  coolant  to 
a  critical  heat  flux  condition.  A  review  of  literature 
concerning  critical  heat  flux  for  water  in  flow  boiling 
indicates  that  many  attempts  have  been  made  to  correlate 
CHP  data  under  diverse  conditions.  Unfortunately,  ex- * 
trapolation  of  any  correlation  beyond  its  data  range  is 
quite  tentative,  as  the  predictions  of  different  corre¬ 
lations  may  vary  widely.  Pertinent  correlations  and 
their  data  ranges  47 pear  in  Table  II-l.  British  Engineer¬ 
ing  Units  will  be  u*  throughout  this  chapter,  since 
most  correlations  are  formulated  using  them. 

To  evaluate  tha  validity  of  CHF  correlations  for  a 
2 

1  kw/cm  heat  flux,  the  constraints  on  mass  flux,  G,  and 


inlet  temperature,  T.  ,  were  defined  for  base  case  con- 

in 

ditions,  and  plotted  in  Figure  II-l.  Reasonable  CHF  pre¬ 
dictions  should  fall  in  a  range  between  the  onset  of  nuc¬ 
leate  boiling  and  the  flow  regime  transition  from  slug  to 

annular  flow.  T.  could  vary  from  freezing  to  saturation 

in 

temperature . 

The  conditions  at  which  the  onset  of  nucleate  boiling 
occurs  just  at  the  end  of  the  channel  may  be  determined 


from  the  equation: 


(1)  G  - 


where; 


4  q"  L  (^) 
If  DC 


(T«t-Tin»  -  +  iTONB 


,  k  (1  1,  '  ^  CD  *4  Aw  -14 

(2J  h  =  .023  0~f)  (-g — i)  (— -)  ,  (Ref. 16), 
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.435  p 


.0234 


|  15.60  p . J 

is  thl  Rchsenow-  Bergles  correlation  for  incipient  nuclea- 
tion  in  British  Engineering  Units  (Ref.  17) . 

Mass  flux  must  be  large  enough  to  avoid  transition 
from  slug  to  annular  flow  in  the  channel,  which  occurs 
at  approximately  an  equilibrium  quality  of  0.1  at  the  base 
case  pressure.  From  the  energy  balance,  a  lower  bound  on 
G  is  found  by  requiring  that  the  equilibrium  quality  at 
the  channel  exit  be  less  than  10%: 


G  ^ 


4  q"  L  (^) 


ir  dFc  (t  -t.  )+.i  ht  ] 
L  P  sat  in  fgj 


The  condition  of  exit  saturation  is  determined  similarly 
and  is  also  plotted  on  Figure  II-l. 

The  predictions  of  various  correlations  were  plotted 
for  the  base  case  as  in  Figure  II-2.  Neither  of  the 
established  physical  mechanism  boundaries  on  mass  flux 
appear  to  come  close  to  predicting  CHF.  A  better  theoreti¬ 
cal  approach  would  be  determination  of  the  point  of  net 
vapor  generation,  following  the  procedure  of  Zuber  and 
Saha (Ref.  54).  This  procedure  stipulates  the  conditions 
under  which  bubbles  first  detach  from  the  channel  walls. 

As  the  mechanism  for  CHF  at  very  high  heat  fluxes  is  pro¬ 
posed  to  be  insulation  of  the  surface  by  a  bubble  boundary 
layer  (Ref.  18)  this  net  vapor  generation  point  is  thought 


i:c  ue  a  good  prediction  of  the  point.  A  p. 


Gl  cue 


Saha-Zubar  equation  appears  on  Figure  II-2,  showing  good 
agreement  with  the  Bernath  and  Gunther  correlations  at 
intermediate  subcoolings.  The  agreement  is  not  as  close 
when  tube  length  is  varied,  but  the  Saha-Zuber  prodedure 
provides  a  CHF  prediction  based  on  a  plausible  theoretical 
mechanism.  It  should  be  evaluated  over  a  range  of  variables 
and  checked  for  agreement  with  available  CHF  data. 

A  comparison  of  the  plots  in  Figure  II-2  with 

corresponding  data  ranges  of  Table  II-l  indicates  that 

the  Bernath  correlation  is  in  general  most  appropriate 

for  the  base  case  conditions.  It  purports  to  be  valid 

over  such  a  wide  range,  however,  that  its  accuracy  is 

only  to  within  30%  (Ref.  18).  Further,  its  range  does 

not  include  particularly  high  subcooling,  where  its  plotted 

behavior  becomes  erratic.  Since  high  subcooling  appears 

to  permit  the  most  efficient  operation,  a  correlation 

suited  to  that  condition  is  desirable.  The  correlation 

of  Rousar  (Ref.  19)  accounts  for  high  subcooling  at 

elevated  pressures  for  tube  lengths  up  to  15  cm.  It  is 

based  on  uniform  circumferential  heat  fluxes  from  1.6 
2 

to  10.3  kw/cm  ,  and  suggests  that  for  very  high  heat 


■c  .  *  - 


,» ••  .  •* 


•fluxes,  local  conditions  govern  CHF,  independent  of 

pressure,  length,  and  diameter.  Work  by  Lowdermilk,  et 
? 

al.  (Ref. 20)  with  thin  tubes  up  to  114  cm  in  length  offers 

2 

data  in  the  1  kw/cm  range,  and  indicates  a  length  and 

diameter  effect.  It  does  not,  however,  include  the 

effect  of  varied  subcooling.  In  order  to  permit  variation 

of  length,  diameter,  and  subcooling  in  parametric  analyses 

to  follow,  a  combination  of  the  Rousar  and  Lowdermilk 

2 

correlations  is  proposed  for  the  1  kw/cm  heat  flux  range. 
The  Rousar  treatment  of  CHF  as  a  local  condition  propor¬ 
tional  to  subcooling  is  applied  to  the  Lowdermilk  correla¬ 


tion,  giving: 


1400  I”  G  AT  .  1  *5 
sub 

,-°V15L  90  _  ^ 


For  a  uniform  axial  heat  flux  the  critical  condition  would 


occur  at  the  tube  exit,  where: 


=  T  ,  -  T .  - 

sat  in 


4  q"  L 


ff  D  c  G 

D 


.  (6) 


Substituting  into  the  above  expression  and  rearranging  gives: 


’crit 


D'°V15  q"'  2 

1400 


~4  q" 
ir 


(T  .  -  T .  ) 
sat  in 


This  expression  closely  corresponds  with  the  Rousar 
correlation  for  the  base  case,  as  shown  in  Figure  II-3. 

The  effects  of  varying  pressure,  length,  and  heat  flux 
on  the  behavior  of  the  Bernath,  Rousar,  and  modified 
Lowdermilk  correlations  are  plotted  in  Figures  II-3, 
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II-4,  and  II-5.  The  modified  Lowdermilk  correlation  is 

seen  i^o  give  more  conservative  estimates  of  CHF  than 

either  Rousar  or  Bernath  at  longer  channel  lengths.  This 

conservatism  will  be  incorporated  into  the  design  window 

to  compensate  for  the  lack  of  specific  data  over  the 

ranges  of  interest.  Figure  II-5  indicates  that  the 

modified  Lowdermilk  expression  is  not  conservative  for 

higher  heat  fluxes.  Design  window  development  at  heat 

2 

fluxes  in  the  range  of  Rousar1 s  data  (q” >  1.6  kw/cm  ) 
should  use  Rousar' s  correlation,  but  with  the  caution  that 
extending  the  channel  length  may  make  the  correlation 
unconservative . 

In  order  to  account  for  experimental  error  and 
component  tolerances,  the  common  design  practice  in  fission 
reactors  has  been  to  establish  a  minimum  critical  heat  flux 
ratio  below  which  the  coolant  should  not  go  during  normal 
operation.  Until  recently,  this  design  CHFR  has  been  1.3, 
meaning  that  the  predicted  CHF  must  always  be  at  least 
1.3  times  larger  than  the  actual  local  heat  flux.  Exten¬ 
sive  experience  with  fission  reactors  has  recently  led  to 
the  lowering  of  this  CHFR  in  fission  reactor  design.  Due 
to  the  uncertainties  inherent  in  the  fusion  program,  the 
safety  factor  of  1.3  will  be  retained  for  this  analysis. 

The  effect  'of  lowering  the  safety  factor  is  examined  in 
Chapter  VI. 


One  means  of  acnieving  higher  heat  loads  is  simply 
9 

increasing  the  coolant  mass  flow  rate.  This  option  is 

constrained,  however,  by  the  unavoidable  pressure  loss 

experienced  during  passage  along  the  coolant  channel. 

Extremely  high  heat  loads  have  been  handled  successfully, 

where  the  heated  length  was  on  the  order  of  1-2  cm,  by 

pumping  coolant  at  very  high  velocity  and  with  large 

pressure  losses  (refs.  21,  22).  This  method  is  employed 

at  NASA- Ames  Research  Laboratory,  where  a  540  psi  pressure 

drop  is  experienced  over  a  10  cm  cooling  length  to  accomo- 

2 

date  heat  loads  as  high  as  12  kw/cm  (Ref.  9).  For  the 
divertor,  somewhat  longer  heated  lengths  may  be  desirable, 
accomodating  lower  heat  fluxes,  requiring  more  reasonable 
pressure  losses,  and  lower  pumping  power  requirements. 

For  the  case  of  single  phase  flow  the  pressure  drop 
to  be  expected  in  the  channel  may  be  approximated  using  the 
Darcy-Weisbach  friction  factor  formula,  taking  the  friction 
factor  to  be  an  empirical  approximation  to  the  Moody  chart 
for  smooth  pipes  at  high  Reynolds  numbers.  Complications 
arise  due  to  the  variation  of  viscosity  with  the  non- 
uniform  temperature  distribution  in  the  flow  cross  section. 
For  uniform  circumferential  heating,  a  viscosity  ratio 
factor,  {U  >  is  commonly  included  in  the  friction 

factor.  Suggested  values  for  the  exponent  v  range  in  the 
literature  from  .  i-i  u^ef .  16.'  „o  .35  (Ref'.  28).  For  single- 


sided  heating,  the  film  temperature  drop  would  vary  around 
the  channel  circumference.  Kays  (Ref.  56)  provides  a  pro¬ 
cedure  for  determining  the  local  heat  transfer  coefficient 
inside  a  tube  under  single-sided  heating.  By  analogy,  an 
expression  for  local  friction  factor  may  be  obtained,  and 
the  average  tube  friction  factor  determined.  For  simplicity, 
a  viscosity  ratio  factor  of  is  employed  here, 

where  U  is  evaluated  at  a  wall  temperature  averaged 
around  the  tube.  Properties  for  water  are  taken  from  Ref. 

23.  The  single  phase  pressure  drop  equations  are  then: 


Ap 


f  L  G 


2  D 


/  gc 


where ; 
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.184  (^b)  '2  fM w  )  ' 

(G  D)  *2  *  Mb  ’ 


,  3xl04  <  Re  ^107,  (9) 


L  =  channel  lenght 
D  =  channel  diameter 
G  =  coolant  mass  flux 
P  =  coolant  density 
yi(b  =  bulk  coolant  viscosity 

JA  -  coolant  viscosity  at  average  wall  temperature 

g  =  units  conversion 
c 

At  high  heat  loads  and  moderate  mass  flux,  subcooled 


boiling  along  the  heated  channel  walls  is  very  possible. 
Several  investigators  provide  some  indication  of  the 


pressure  losses  to  be  expected,  among  them;  Owens  and 
Schrock  (Ref.  24),  Tarasova  (Ref.  25),  Mendler  (Ref.  26), 
Martinelli  and  Nelson  (Ref.  27),  and  Dormer  and  Bergles 


(iv-f.  28).  Wide  var.umc  • 


in  pressure  drops  and  test 


conditions  exists,  however,  suggesting  more  design  specific 
experimentation  before  taking  full  advantage  of  the  sub¬ 
cooled}  boiling  condition. 


The  study  of  first  wall  cooling  by  M.A.  Hoffman,  et 
al.  (Ref.  12),  concluded  that  the  Mendler-Martinelli- 
Nelson  model  was  most  appropriate  for  their  design  case, 
resulting  in  a  tripling  of  the  effective  friction  factor 
over  the  nonboiling  factor  near  the  end  of  their  tube. 

The  Owens-Schrock  model  predicted  a  thirteen-fold  in¬ 
crease  in  friction  factor  for  the  same  subcooled  boiling 
conditions.  The  experimental  work  of  Dormer  and  Bergles 
(Figure  II-6)  records  the  behavior  of  total  pressure  drop 
for  various  tube  L/D  ratios  as  the  subcooled  boiling 
regime  is  entered.  For  L/D  values  of  50  and  100,  total 
tube  pressure  drops  in  subcoolev  boiling  more  than  six 
times  greater  than  the  corresponding  adiabatic  flow 
pressure  drop  were  recorded.  For  shorter  tubes  (L/D  =  25) , 
pressure  drops  had  reached  only  about  3  times  their  adia¬ 
batic  value  when  onset  of  critical  heat  flux  occurred. 

For  tubes  with  L/D  =  10,  the  higher  wall  shear  common  to 
the  entrance  region  may  be  expected  to  increase  the  total 
pressure  drop  by  up  to  50%  (Ref.  16),  so  that  total  pres¬ 
sure  drops  up  to  4  to  5  times  larger  than  the  non-boiling 
AP  may  be  expected  when  fully  developed  subcooled  boiling 
occurs . 

suhcooled  boiling  pressure  drop  in 


Evaluation  of  tho 


the  10  cm  long  base  case  tube,  using  the  Mendler  approach 
(Ref.  26} ,  indicates  that  only  at  very  low  local  subcooling 
does  €he  subcooled  boiling  pressure  drop  exceed  the  pre¬ 
dicted  isothermal  pressure  drop.  As  is  evident  in  Figure 
II-6,  the  favorable  viscosity  effect  overcomes  the  boiling 
effect  shortly  after  onset  of  boiling.  Only  when  the  heat 
flux  on  the  tube  becomes  a  significant  fraction  of  that 
which  would  cause  the  exit  to  be  saturated  (under  the  same 
flov  conditions)  does  the  subcooled  boiling  effect  greatly 
raise  the  pressure  drop.  For  longer  tubes,  both  the  Dormer 
data  and  the  Mendler  formula  predict  a  more  rapid  increase 
in  pressure  drop  upon  onset  of  boiling,  necessitating  more 
careful  analysis. 

From  Figure  II-4  it  is  apparent  that  at  sufficiently 

high  subcoolings,  adequate  cooling  can  be  provided  without 

2 

boiling  for  a  1  kw/cm  heat  flux,  in  tubes  up  to  1  meter 

7  2 

long,  with  a  mass  flux  on  the  order  of  10  lbm/hr-ft 
The  non-boiling  pressure  drop  over  a  1  meter  long  tube  at 
this  mass  flux  is  less  than  3  psi,  so  that  even  a  six¬ 
fold  increase  in  Ap  due  to  subcooled  boiling  should  not 
be  a  major  problem.  Only  at  lower  inlet  pressures,  higher 
heat  loads,  and  longer  lengths  would  the  subcooled  boiling 
pressure  drop  become  critical. 

The  ratio  of  pumping  power  required  for  coolant 
circulation  to  thermal  power  extracted  by  the  coolant  is 


a  frequently  used  measure  of  cooling  system  efficiency. 

5  Pumping  Power  m  ap//> 

W  =  -  =  - 

p  Thermal  Power  me  AT,  - 

p  1-2 

AP 

=  -  (10) 

P  °P  4TI-2 

Based  on  overall  power  conversion  efficiency  in  a  fusion 
reactor  blanket,  Fraas  (Ref.  29)  suggested  a  design  limit 
for  this  ratio  of  2%.  Since  the  divertor  would  handle 
only  a  small  fraction  of  the  total  plasma  power,  it  could 
be  run  at  a  much  higher  ratio,  without  very  seriously  im¬ 
pacting  overall  plant  efficiency.  However,  since  the 
occurrence  of  subcooled  boiling  may  increase  the  pressure 
drop  up  to  six  times,  the  2%  limit  will  be  used  for  the 
design  window,  with  Ap  calculated  for  non-boiling  flow. 

W  =  — — P— °-n-  <  .02  .  (ID 

P  J>  c  AT  _ 

*  p  1-2 

Figure  I I- 7  indicates  the  curve  for  pumping  power  ratio 
=  2%  in  the  nonboiling  and  boiling  regions.  The  small  area 
included  below  the  extension  of  the  nonboiling  W^=  2%  into 
the  subcooled  boiling  region  represents  a  "flagged"  region 
in  the  design  window.  Operating  points  here  require 
special  analysis  for  subcooled  boiling  pressure  drop,  using 
the  Mendler  approach. 

D.  Heat  Transfer  Coefficient 

•;,ie  susccptiL  ..lit/  cf  t;.  t  ...i  trials  tc  -  nach 

temperatures  at  which  material  integrity  diminishes  de- 


mands  that  good  heat  transfer  from  channel  wall  to  fluid 
be  achieved.  Subcooled  boiling  and  other  mechanical 

methods  enhance  fluid  mixing,  thereby  increasing  the  heat 
transfer  coefficient.  For  the  design  window  a  single  phase 
heat  transfer  coefficient  by  McAdams  (Ref.  19)  will  be  used 
to  develop  the  fluid-wall  interface  temperature  profile. 

Nu  =  .023  (Re)  ‘  8  (Pr)  ‘  4  Ub/«w)  * 14  (12) 

The  occurrence  of  subcooled  boiling  would  make  this  equa¬ 
tion  conservative.  Further,  the  steep  temperature 
gradient  at  the  wall  in  the  entrance  region  would  enhance 
beat  transfer  there,  further  increasing  the  conservatism 

for  short  tubes.  For  base  case  conditions  and  mass  flux 
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on  the  order  of  10  lbm/hr-ft  ,  this  equation  predicts  a 
film  temperature  rise  of  "210  °F.  Even  with  low  inlet 
temperature  and  subcooled  boiling,  therefore,  the  inner 
wall  temperature  will  be  over  300  °F,  posing  a  high  temp¬ 
erature  materials  challenge. 

E.  Swirl  and  Mixed  Flow  Schemes 

Extensive  work  has  been  documented  on  the  advantages 
realized  in  coolant  channel  performance  when  internal  de¬ 
vices  are  employed  to  induce  swirl  flow.  The  effect  is  to 
strip  away  nucleating  bubbles  and  force  denser  (colder) 
coolant  to  the  outside  wall  where  heat  transfer  is  occur¬ 
ring.  Correlations  for  critical  heat  flux,  pressure  drop, 
and  heat  transfer  coefficients  in  such  cases  are  presentee 


by  Garabill,  et  al.  (Ref.  7).  Approximately  1.5  to  1.8 

times  higher  CHF  can  be  realized  over  straight  flow  at  the 
c 

expense  of  1.4  to  1.8  times  higher  pressure  drop.  The 
heat  transfer  coefficient  is  improved  over  straight  flow 
by  up  to  85%  for  constant  flow  rate  (Ref.  30).  In  terms 
of  available  design  window  space,  the  benefits  of  higher 
CHF  are  roughly  offset  by  the  penalty  of  higher  pumping 
power.  Swirl  flow  appears  attractive  in  easing  material 
temperature  restrictions  through  increased  film  heat 
transfer  coefficient.  Unless  operation  near  a  material 
temperature  limit  is  required,  the  use  of  swirl  flow  does 
not  appear  necessary  for  the  base  case  conditions.  From 
a  particle  handling  viewpoint,  swirl  flow  offers  a  means 
of  controlling  diffusion  rates  through  channel  wall  temper¬ 
ature  control. 

F.  Effect  of  Asymmetric  Heating 

Collier  (Ref.  31)  summarizes  data  regarding  the  effects 
of  asymmetric  circumferential  heating.  For  highly  sub¬ 
cooled  boiling,  he  recommends  treating  the  occurrence  of 
burnout  as  a  localized  condition.  This  implies  that  burn¬ 
out  occurs  where  the  incident  flux  is  normal  to  the  channel 
wall,  at  the  same  heat  flux  as  would  cause  burnout  for 
the  uniform  heating  case.  This  is  clearly  conservative, 


as  no  allowance  is  made  for  fluid  convection,  conduction 
around  the  tube  circumference,  or  conduction  in  the  fluid. 


As  a  conservative  limit,  the  uniformly  heated  CHF  value 
will  be  considered  the  design  constraint. 

The  case  of  swirl  flow  virtually  negates  the  effect 
of  asymmetry,  since  the  entire  fluid  column  is  rotated 
past  the  hot  spot.  Kim,  et  al.  (Ref.  32)  predict  that  CHF 
values  twice  as  high  as  for  uniform  heating  may  be  achieved 
utilizing  a  swirl  flow  for  cooling  a  one-side-only  heating 
situation. 
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III.  SURFACE  INTERACTIONS 


A.  Spluttering  and  Plasma  Contamination 

The  problem  of  impurity  sputtering,  which  the  divertor 
is  intended  to  reduce  at  the  tokamak  first  wall,  is 


particularly  acute  at  the  divertor  target  itself.  The 
ability  of  energetic  ions  to  erode  the  solid  target  surface 
upon  impact  is  a  major  limit  on  the  target  design.  The 
rate  at  which  the  material  will  erode  may  be  expressed  to 


first  order  by  the  formula: 


L  t 

Tr 


2 

i 


S.  J. 
1  1 


C  ( 


M 


N. 


(13) 


where;  At  =  change  in  thickness 

tT  =  change  in  time 

=  sputtering  coefficient  (atoms/ion) 

J.  =  ion  flux  at  divertor  surface 

l 

C  =  capacity  factor  of  target  (%  of  year 
exposed  to  target) 

M  =  target  material  atomic  weight 

P  =  target  material  density 

Na  =  Avogadro ' s  Number 

This  relation  may  be  used  in  comparing  candidate 


materials  for  use  as  target  materials.  The  distribution 


of  particle  fluxes  and  energies  across  the  target  surface 


must  be  considered  to  account  for  accelerated  erosion  in 


regions  of  peak  flux.  Lacking  such  information,  this  work 
will  apply  the  average  particle  fluxes  to  impose  a  com- 

rr-  -ative  conr^rair*-  mater?  .*  “  th cVness  and  1-  fetir 

Sputtering  coefficients  are  taken  from  Refs.  33  and  34. 


The  energy,  mass,  and  angle  of  incidence  of  impinging 
ions  have  significant  effect  on  the  local  sputtering  coef- 

t 

ficient.  Figure  III-l  shows  the  typical  dependance  of 
light  ion  sputtering  coefficient  on  incident  energy  and 
mass.  The  peak  of  each  curve  occurs  in  the  range  1-2  keV. 

The  sputtering  coef¬ 
ficient  varies  with  incident 
angle  according  to: 

si(0)  =  '  (Ref*  52) 

where  8=0  constitutes  normal 

incidence. 

Since  a  negative  sheath 
potential  will  exist  at  the 
surface  of  an  operating  di¬ 
vertor,  ions  will  be  accel¬ 
erated  abruptly  toward  the 
target.  This  acceleration 
should  raise  the  ion  energy 
from  the  plasma  edge  temperature  of  .3  to  .5  keV  up  to  a- 
bout  1.3  keV  (Ref.  53).  Further,  it  should  make  the  angle 
of  incidence  nearly  vertical.  Normal  incidence  at  1.3  keV 
is  assumed  for  this  study.  The  effect  of  improving  the 
sputtering  climate  is  examined  in  Chapter  VI.  It  should  be 
noted  that  helium  and  heavier  impurities  have  not  been  in¬ 
cluded  as  sputtering  agents  for  this  analysis. 

T.-.ble  III-l  su:.;.wi.:  :s  the  key  sputtering  variables 
and  their  product  for  the  eight  candidate  materials. 


Figure  III-l:  Typical 
Sputtering  Coefficient  Data 


As  shown  xn  figures  III- 2  and  III- 3,  the  sputtering  life¬ 
time  of  a  given  target  material  increases  linearly  with 

9 

its  thickness.  This  suggests  a  minimum  thickness  per¬ 
missible  for  a  specified  target  life. 

Considerable  investigation  is  warranted  to  determine 
to  what  degree  any  redeposition  of  surface  material  may 
be  expected  on  the  exposed  target  surface.  It  would  seem 
implausible  for  such  assistance  to  occur  in  regions  of  peak 
flux,  though  its  demonstration  would  certainly  be  welcome. 
None  is  assumed  here.  The  effect  of  self  sputtering  must 
also  be  addressed;  a  destructive  cascade  being  possible 
for  the  heavier  materials. 

Mere  comparison  of  sputtering  rate  is  not  sufficient 
measure  of  a  candidate's  usefulness  as  a  target  material. 

The  destiny  and  effect  on  the  fusion  plasma  of  the  sput¬ 
tered  impurity  must  also  be  considered.  A  limit  on  the 
impurity  concentration  in  the  plasma,  above  which  ignition 
is  impossible,  is  provided  by  the  expression  (Ref.  35) : 

Imax  -  50  *  W1  Z  <  <14> 

where;  I  =  maximum  allowable  impurity  (  %  ) 

max  r  J 

Z  =  target  Z  number 

Assuming  a  plasma  volume  characteristic  of  next  gen- 

.  3 

eration  reactors  (  400  m  )  allows  determination  of  a 
maximum  allowable  impurity  concentration.  Further, 
assuming  a  divertor  duct  geometry  and  pumping  system  capable 
of  trapping  99%  of  the  target  generated  impurities  estab- 


lishes  an  estimate  of  the  duration  of  target  operation 
permissible  before  impurities  from  the  target  penetrate 
the  plasma  core  in  sufficient  quantity  to  extinguish  the 
burn.  These  results  are  tabulated  in  Table  III-2  for 
nine  materials.  Such  results  are  tentative,  since  the 
divertor  pumping  efficiency  and  the  tendency  for  impurities 
at  the  plasma  edge  to  migrate  into  the  plasma  core  are  at 
present  uncertain.  The  comparison  of  Table  III-2  is 
valuable,  however,  for  comparing  material  alternatives. 

A  relative  figure  of  merit  may  be  devised  by  simply 
dividing  the  sputtering  rate  by  the  allowable  burn  time. 

The  ranking  in  Table  III-3  results,  indicating  the 
undesirability  of  copper  and  tungsten,  and  the  advantages 
of  beryllium  and  graphite  as  plasma  interface  materials. 

B.  Internal  Erosion 

The  thickness  of  the  target  material  may  also  be  dim¬ 
inished  by  erosion  and  corrosion  on  the  coolant  side. 
Corrosion  is  a  relatively  slow  process  for  most  materials 
under  study,  and  should  have  little  effect  relative  to  the 
external  sputtering  rate.  Common  corrosion  rates  are  at 
least  three  orders  of  magnitude  smaller  than  calculated 
sputtering  erosion  rates  (Ref.  36).  Mechanical  erosion 
by  the  rapidly  flowing  coolant  may,  however,  be  a  concern. 
Experience  with  high  mass  flow  boilers  and  studies  of  rain 
damage  to  high  Performance  aircraft  have  indicated  that 
liquid  water  travelling  at  speeds  in  excess  of  50  meters 
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per  second  will  cause  mechanical  damage  to  metal  surfaces. 

Coolant  velocities  of  less  than  30  meters  per  second  should 

*  2 

be  adequate  for  a  1  kw/cm  heat  flux,  so  direct  erosion 
should  not  occur.  The  incidence  of  surface  boiling  may 
generate  regions  of  bubble  collapse  with  erosion  rates  com¬ 
parable  to  those  experienced  in  cavitation  erosion.  Such 
erosion  has  a  velocity  below  which  erosion  is  undetectable. 
Above  this  threshold,  the  erosion  rate  increases  approxi¬ 
mately  as  the  velocity  to  the  sixth  power  (Ref.  37) . 

Such  erosion  would  place  a  distinct  limit  on  the  flow 
velocity  and  hence  the  allowable  heat  load. 

The  main  damage  mechanism  in  cavitation  erosion  occurs 
when  voids  in  the  fluid  suddenly  collapse,  generating  both 
intense  shock  waves  and  liquid  microjets  capable  of 
pitting  or  work  hardening  nearby  surfaces  until  fatigue 
failure  and  material  loss  occur.  Voids  generated  by  boil¬ 
ing  on  a  hot  surface  cannot  collapse  at  the  surface  in  the 
same  way  as  cavitation  bubbles,  since  vapor  is  always 
being  generated  by  the  addition  of  heat.  In  subcooled 
boiling,  the  bubble  collapse  should  generally  occur  away 
from  the  heated  walls  and  at  a  slower  rate  than  in  cavita¬ 
tion.  Thus  any  internal  erosion  should  be  downstream  from 
the  heated  surface,  and  of  lesser  intensity  than  that  due 
to  cavitation.  Particular  attention  will  therefore  be 
necessary  where  the  flow  direction  changes  downstream  of 
the  heated  section  to  insure  that  failure  by  internal 


erosion  does  not  occur. 

Experiments  with  cavitation  erosion  on  aluminum 
surfafees  indicated  that  no  appreciable  damage  was  sus¬ 
tained  at  water  velocities  less  than  12  m/sec  (Ref.  37). 
Since  adequate  cooling  can  be  provided  at  high  subcooling 

with  velocities  less  than  10  m/sec  (Fig  II-7) ,  internal 

2 

erosion  should  not  be  a  problem  for  a  1  kw/cm  heat  flux 
case.  Should  higher  velocities  be  employed  to  offer 
protection  against  plasma  disruptions,  the  possibility  of 
cavitation  damage  must  be  specifically  examined.  TZM  is 
expected  to  have  particularly  good  resistance  to 
cavitation  damage  in  water,  since  its  oxide  is  one  of  very 
few  known  to  resist  wear  better  than  its  base  metal  (Ref. 


Be 

4 

.036 

9.012 

1.848 

.176 

C 

gr 

6 

.019 

12.011 

2.10 

.109 

A1 

13 

.043 

26.982 

2.70 

.429 

Ti 

22 

.011 

47.90 

4.54 

.116 

V 

23 

.014 

50.94 

6.11 

.117 

Cu 

29 

.101 

63.546 

8.96 

.716 

Nb 

41 

.0046 

92.91 

8.57 

.050 

Mo 

42 

.0066 

95.94 

10.22 

.062 

Table  III-l:  Key  Sputtering  Characteristics  for  Eight 
Candidate  Elements 
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a.  Max  impurity  concentration  for  ignition: 

I  =  50  x  10-'1  Z  (Ref.  35) 
max 

b.  ETF  plasma:  6.4  x  10^2  particles. 

c.  Sputtering  rate  =  J.^. 

d.  Assume  1%  of  sputtered  material  gets  into  plasma. 

(99%  is  pumped  or  condensed.) 

e.  Plasma  burn  time  limit  due  to  impurity  contamination 

max  impurity  allowed 
sputtering  rate  @  1  kw/cm^ 


f.  Results 

W  -  .0082  sec  »  Q  A1  -  234  sec 

Mo  -  1.94  sec  TJ  c  Be  -  2178  sec 

Cu  -  2.46  sec  C  -  2610  sec 

Ti  -117.  sec  gr 


Table  III-2:  Effects  on  Plasma  Burn  Times  of  Sputtered 
Impurities  from  Candidate  Target  Materials. 


Figure  of  Merit 


sputter  rate  ,  mm  . 
burn  time  yr-sec'7 


W  -  963.4  Nb 
Mo  -  254.5  A1 
Cu  -  27.8  V 


10.6 

Ti  - 

.838 

1.58 

Be  - 

.071 

1.41 

Cgr” 

.037 

Table  III-3:  Comparative  Figure  of  Merit  for  Candidate 
Materials  Due  to  Sputtering  and  Plasma 
Contamination  (Small  Numbers  Preferred). 


IV.  THEPMO-MECHANICS 


A.  Pressure  Stress 

The  containment  of  internal  coolant  pressure  imposes 
a  stress  constraint  which  demands  that  channel  wall  thick¬ 
ness  always  exceed  a  certain  minimum.  For  cases  in  which 
the  wall  may  be  considered  the  membrane  of  a  thin  walled 
cylinder,  the  stress  intensity  is  easily  determined  as 
the  hoop  stress  minus  the  radial  stress: 

S  =  (  p  r  /  t  )  -  (  -p  ) .  (15) 

The  greatest  allowable  stress  intensity  within  a 

material  is  established  by  the  ASME  Code  as  S  .  ,  and  must 

mt 

be  less  than  1/3  of  the  ultimate  strength  and  less  than  2/3 
of  the  yield  strength.  The  minimum  wall  thickness  allowed 


For  thick  walled  pressure  vessels,  the  stress  at  any 


point,  n,  in  the  cylindrical  cross  section  is  given  by: 


(17) 


For  internal  pressurization,  the  tangential  stress  is 
tensile,  the  radial  stress  compressive,  and  the  axial 


which -is  maximum  at  the  inside  surface  of  the  cylinder. 


n=a.  This  requires  that: 

2 


2  b 


b2  -  a 


2  (  P  )  *  Smt  . 


(19) 


Substituting  a=r,  b=r+t  leads  to  determination  of  the 
minimum  thickness  for  this  case: 

1 

-  )  r  .  (20) 


t  .  =  (  -1  + 

mm 


yi  -  2 p/s 


mt. 


This  is  slightly  larger  than  the  minimum  thickness  with 
the  thin  wall  assumption  and  will  be  used  when  f>0.1  r  . 


B.  Thermal  Stress 

Heat  loads  of  the  magnitude  expected  at  the  divertor 
target  introduce  steep  temperature  gradients  in  target 
material.  Cyclic  operation  induces  large  temperature 
changes  between  plasma  burn  and  chamber  rejuvenation 
periods.  The  result  of  these  conditions  is  sizeable 
thermal  stresses  which  dominate  the  total  stress  calcula-' 
tions. 

Careful  design  is  necessary  to  minimize  the  external 
constraints  on  material  expansion  under  heat  load.  Self¬ 
constraint  is,  however,  unavoidable  and  must  be  considered 
in  determining  the  peak  stress  to  which  the  material  is 
subjected.  For  sections  of  a  plate  or  thin  walled  cylin- 
d~  ical  tu..5jt  sufficiently  far  from  the  edge  attachments, 
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under  s*'  •  ■■  v  state  condition-:,  the  temperature  profile  may 
be  assumed  to  be  nearly  linear  from  the  high  temperature 
heate<3  surface  to  the  low  temperature  cooled  surface.  For 
a  plate  constrained  against  bowing  and  for  a  cylindrical 


tube,  the  peak  thermal  stress  will  be: 


peak 


E  °<  AT 

2  (  l-tf) 


(Ref.  40) 


compressive  on  the  hot  surface,  tensile  on  the  cold.  The 
difference  in  temperature  may  be  expressed  in  terms  of  the 
heat  flux,  wall  thickness  and  conductivity: 

AT  =  .  ,72) 


Thus  a  rimple  expression  relating  the  maximum  thermal 
stress,  the  heat  flux,  and  wall  thickness  is: 

E  a"  t  ^ 

c rth  =  - : •  (23) 

peak  2  k  (  1-  ))  ) 

Thermal  stresses  at  edge  connections  become  difficult  to 
determine  analytically  and  require  specific  study  for 
particular  designs.  Except  in  the  cases  of  a  free  edge 
or  specific  designs  in  which  the  edge  attachments  deflect 
readily,  edge  connections  will  cause  a  rise  in  stress  over 
the  "infinite  wall"  case.  Equation  23  is  therefore  an 
optimistic  estimate  of  the  thermal  stresses  to  be  encount¬ 
ered. 


C.  Fatigue 

The  L«  ‘.e  case  heat  load  cycle  considered  in  this 


•.*  \s  cor 


foil  owe.  1 
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second  quiescent  rejuvenation  period.  As  a  result,  over 
105  cycles  are  likely  during  one  year‘s  operation,  intro- 
ducing  the  likelihood  of  fatigue  failure  as  a  controlling 
design  consideration.  This  may  be  treated  using  semi- 
empirical  fatigue-life  curves  for  several  candidate  mater¬ 
ials  in  conjunction  with  ASME  Code  rules  for  peak  stress 
analysis.  The  number  of  cycles  anticipated  dictates  that 
the  fatigua  will  be  strain-cycle  controlled. 

Analysis  of  the  time  dependent  temperature  profile 
through  the  target  material  indicates  the  temperature  and 
strain  ranges  through  which  the  material  is  cycled. 
Considering  a  semi-infinite  body  whose  surface  is  exposed 
to  periodic  temperature  cycling,  the  fraction  of  the  surface 
temperature  fluctuation  experienced  at  a  depth,  x,  is 

Vo  -  «*p  <  -  >  -  <2«> 

where  (J  is  the  period  of  oscillation,  and  afc  is  material 

thermal  diffusivity.  This  suggests  that  thermal  strains 

inside  the  material  might  be  less  than  calculated  using 

steady  state  equations,  due  to  the  oscillating  heat  load. 

The  divertor  target  is  not  infinite,  however,  and  is 

more  closely  modelled  by  a  slab  with  fixed  temperature,  T  , 

c 

on  one  side,  and  a  surface  heat  flux  on  the  opposite  side, 

q"  ,  which  follows  a  periodic  step  function  as  in  Figure 
s 


IV- 1. 


99  KS 


IX  to 


Time  (yet) 

Figure  IV-1 


The  time  required  to  reach  an  essentially  steady  state 
t.iv.'Tpexaturc  profile  after  a  step  change  may  be  determined 
,'J.:.om  the  integral  technique  as  (Ref.  41)  : 


For  practical  slab  thicknesses  (1  to  10  mm)  and  the  eight 
materials  considered  in  the  analysis,  these  response  times 
are  at  least  an  order  of  magnitude  less  than  fifteen 
seconds.  Fatigue  strain  ranges  are  therefore  calculated 
based  on  full  temperature  cycling  between  a  steady  state 
burn  condition  and  a  steady  state  off  condition. 

Data  for  strain-cycle  fatigue  life  curves  is  found 
in  References  42,  43,  and  44.  Peak  stresses  determined  as 
a  combination  of  the  pressure  stress  and  thermal  stress  on 
the  structure  are  converted  to  equivalent  strains  by  use 
of  the  material  modulus  of  elasticity,  E.  The  cyclic 
component  of  these  peak  strains  must  be  kept  below  &<max» 
that  strain  which  will  cause  fatigue  fai lure  over  the  life¬ 
time  snecified.  Sub**  i.tutin<-'  -  ,  =  E  (  irto 

.  -  vi  .  . 

equation  23  and  rearranging  gives  a  limit  on  wall  thickness 


'W~ 
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such  that: 


* 


t  * 


2  k  (  1-  V  )  A€mav 
_ max 

rV  q" 


(26) 


This  is  in  conflict  with  the  thicker  channel  walls  favored 
by  the  sputtering  consideration,  and  leads  to  an  optimum 
thickness  for  either  (1)  maximum  allowable  heat  flux  for 
a  specified  lifetime,  or  (2)  maximum  life  at  a  specified 
heat  load.  Determination  of  these  optimum  thicknesses  is 
undertaken  in  Chapter  V,  and  is  the  major  factor  in  the 
design  limitations  reflected  in  the  conclusions  of  Chapter 


VII. 


D.  Thermal  Shock 


During  the  transient  periods  of  burn  initiation  and 
quenching,  a  non-linear  thermal  gradient  will  exist  within 
the  target  material.  This  will  introduce  additional  mech¬ 
anisms  of  fatigue  failure.  The  extremely  steep  temperature 
gradient  at  the  heated  surface  immediately  after  burn 
initiation  may  result  in  cracking  of  the  surface,  with 
possible  spalling  and  accelerated  sputter  damage.  The 

5 

frequently  quoted  thermal  shock  parameter,  crnay  (k  cp  p )  /E«*. 

gives  a  measure  of  material's  relative  survival  in  such  an 
environment  (Ref.  45).  Work  by  Schivell  and  Grove  (Ref.  46) 


indicates  the  greatest  magnitude  of  the  additional  stress 
due  to  thermal  shock  on  a  flat  plate  will  be  -1/5  of  the 
maximum  steady  state  thermal  stress,  and  will  occur  at  a 

wit  ■  rom  the  keae  ;  .  .  ,'uce  approximately  4/1. <  of  the 


total  thickness.  This  location  does  not  coincide  with 
that  of  the  peak  steady  state  thermal  stress  (adjacent  to 
the  coolant) ,  but  the  additive  effects  were  thought  to  be 
significant.  The  thermal  shock  stress  on  the  inside  layer 
is  ~1/10  of  the  steady  state  stress.  Experimentation  under 
cyclic  high  heat  load  conditions  will  be  necessary  to 
determine  actual  effects.  For  this  analysis  the  shock 
stresses  are  assumed  to  consume  a  fraction  of  the  fatigue 
life  equal  to  N/Ng  ^ ,  where  N  is  the  number  of  cycles  en¬ 
dured,  and  ^  is  the  number  of  cycles  to  failure  at  0.1 

times  the  peak  steady  state  thermal  stress.  For  all  the 
materials  considered,  the  shock  stress  fraction  of  fatigue 
life  is  found  to  be  negligibly  small. 


E.  Temperature  Limits 

All  materials  are  subject  to  limitations  in  operating 
temperature,  necessitating  careful  design  to  prevent  over¬ 
heating.  Though  fatigue  life  and  sputtering  are  here  found 
to  govern  design  wall  thicknesses,  the  temperature  restric¬ 
tions  posed  by  materials  characteristics  form  a  secondary 
design  limitation  which  bounds  the  range  of  potential 
improvements . 

The  melting/sublimation  temperature  limit  is  inescap¬ 
able.  In  normal  operation  it  would  occur  first  on  the 
heated  surface  at  the  outlet  end  of  the  coolant  flow,  as¬ 


suming  a  uniforn  axial  heat  flux.  The  limiting  equation  is 


(27 


melt 


“  Tin-  a"’J V+t)  m<l+t/r)  +  q"/h  +  (~^)  2Lq"//TGc^r" , 


1-2' 


the  coolant  temperature  rise. 


where: 

2Lq"  (— — )/  fTGcpr  —  AT 


*  q"/h  =  with  h  as  in  equation  2,  and 

3 — i^tj_  ln(l+t/r)  =  ATwan'  the  temperature  differ- 

through  the  channel  wall  (Ref. 16). 


Most  materials  display  a  distinct  degradation  in 
strength  above  a  critical  temperature,  which  may  pose  a 
design  limit  if  the  material  supports  a  structural  load. 

The  limiting  equation  then  would  be  similar  to  that  above 
with  Tme2t  replaced  by  For  the  divertor  application 

adequate  strength  for  coolant  containment  is  provived  by 
the  inner  cooler  layers  of  material.  Therefore,  degradation 
of  strength  in  the  outer  layers  is  not  considered  limiting. 

Coolant  and  plasma  compatibility  with  structural 
materials  imposes  temperature  limits  as  well.  For  maximum 


coolant  temperature  permissible  to  avoid  excessive  corrosion 
«.  .  (28) 


corr 


T  -T.  >  q"/h  +  (£±i)  2Lq"/VGc  r  =  AT...  +AT. 

corr  in  “  ^  r  p  film  1- 


Known  corrosion  temperature  limits  are  tabulated  in  Figure 
IV- 2. 


F.  Graphical  Materials  Data 

Employing  materials  properties,  equations  26,  27,  and 
28  may  be  plotted  on  heat  flux  vs.  wall  thickness  axes  for 
candidate  materials  as  shown  is  Figures  IV-3a  to  IV-3h. 

Two  fatigue  curves  are  plotted,  representing  different 
■  of  thermal  _  i*.  ■  Lo  failure.  The  pair.;  of  curves 


for  melting  and  property  degradation  temperature  limits 
represent  the  range  of  these  limits  resulting  from  changes 
in  tl$e  heat  transfer  coefficient  due  to  varied  mass  flux 
of  coolant.  Safe  operating  conditions  with  respect  to 
each  material  constraint  lie  below  the  respective  constraint 
line.  To  survive  105  thermal  cycles,  the  wall  thickness 
and  heat  load  must  plot  below  the  10 ^  cycle  fatigue  curve. 

If  fatigue  cracking  and  property  degradation  can  be  accepted, 
(as  in  a  protective  tile  bearing  no  loads) ,  points  below 
the  melt  limit  would  be  acceptable. 


s 

Element 

Gas  Corrosion 

Liquid  Water 
Corrosion 

Liquid  Metal 
Corrosion 

Beryllium 

-Does  not  react  w/ 
hydrogen  up  to 
2050°F 

-Attacked  by  steam 
above  500°F 

-Good  resistance 
to  water  up  to 
500°F 

-Good  resistance 
to  liq.  metals 
-Resists  lithium 
to  1110°F 

Aluminum 

-Hydrogen  doesn't 
diffuse  and  has 
no  corrosive 
effect 

-Good  resistance 
-Protective  film 
forms  in  water 

Titanium 

-Hydrogen  absorbed 
interstitially 
above  600°F 

-Excellent 
resistance  to 
600°F 

-Limited  resist¬ 
ance  to  lithium 
above  147 0°F 

Vanadium 

-Hydrogen  embrittles 
above  570°F 

-Modest 
corrosion  in 
water 

-Poor  resist¬ 
ance  in  general 
-No  tests  on 
lithium 

Copper 

-Hydrogen  em¬ 
brittles  above 

650°F 

-Excellent 
resistance  to 

water 

Niobium 

Easily  oxidized 
-Reacts  easily  w/ 
hydrogen 

-Good  resistance 
to  560°F 

-Good  lithium 
resistance  to 
1650°F 

Molybdenum 

-Oxidizes  rapidly 
above  1000°F 

-Limited 
applications 
in  degassed 
water  to  600°F 

-Good  lithium 
resistance  to 
1830°F 

Graphite 

-Rapid  methane 
production  from 
670-1310°F 
-Rapid  acetylene 
production  above 
2100°F 

• 

Figure  IV-2:  lIaterJ.il  Corrosi  Summary  ORef.  36) 
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TZM  Molybdenum  Alloy 


Figure  IV-3as  Fatigue,  Property  Degradation, 

and  Melting  Limit  Curves  for  TZM 
(Area  Below  Curves  Acceptable) 


(mm) 


Beryllium 


IV-3d:  Fatigue,  Property  Degradation, 

and  Melting  Limit  Curves  for 
Beryllium. 

(Area  Below  Curves  Acceptable) 


FT 
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q"  (  hw/ 'rrT ) 


(mm) 


Sintered  Aluminum  Powder 


Figure  IV- 3f:  Fatigue ,  Property  Degradation, 
and  Melting  Limit  Curves  for 
Sintered  Aluminum  Powder 
(Area  Below  Curves  Acceptable) 


qM  (kw/rn  ) 
* 


Ti-6A1-4V  Titanium  Alloy 


Figure  IV- 3h:  Fatigue,  Property  Degradation, 
and  Melting  Limit  Curves  for 
Titanium  Alloy  (Ti-6A1-4V) 
(Area  Below  Curves  Acceptable) 
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V.  EVALUATION  OF  DESIGN  OPTIONS 

A.  Possible  Geometries 

A  wide  variety  of  solid  target  configurations  are 
possible.  The  range  of  geometries  depicted  in  Figure  V-l 
has  been  considered  in  this  investigation.  Each  may  be 
evaluated  by  applying  minor  modifications  to  the  design 
equations  already  discussed.  For  analysis  purposes,  the 
concepts  may  be  grouped  into  four  more  basic  categories; 
Simple  Tubes  (A,B) ,  Homogeneous  Flat  Plates  (C,D) ,  Compos¬ 
ite  Tubes  (E) ,  and  Armored  Plate  Designs  (F,G,H) . 

B.  Materials 

1.  Functional  Requirements 

Identification  of  suitable  materials  for  divertor 
target  application  is  siruilar  to  the  search  for  tokamak 
first  wall  materials.  Strength  and  ductility  are  desirable 
to  handle  the  thermal  and  pressure  loads  inherent  in  the 
system.  Chemical  compatibility  with  other  materials  in 
the  machine,  particularly  coolant  and  plasma,  must  be  con¬ 
sidered.  Further,  for  materials  exposed  to  the  plasma, 
resistance  to  hydrogen  sputtering  and  self-sputtering  must 
be  sought.  Minimization  of  plasma  quenching  effects  due 
to  impurity  penetration  of  the  plasma  favors  low  Z  materials. 
Resistance  to  thermal  shock  effects,  fatigue  damage,  and 
radiation  damage  are  additional  desirable  characteristics. 

Materials  offering  some  of  the  desirable  characteris- 


.  '*• 
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tics  are  numerous,  and  v.st  potential  improvements  through 
alloying  may  be  available.  An  aggressive  program  of  mat- 
erials  evaluation  is  warranted  and  ongoing  to  determine 
materials  which  meet  all  of  the  many  criteria  imposed  on  the 
divertor  target.  As  a  basis  for  analysis,  eight  structural 
materials  representing  a  range  of  material  characteristics 
are  selected;  Copper (ZAC-2 ) ,  Sintered  Aluminum  Powder, 
Beryllium,  Graphite,  Molybdenum  Alloy (TZM),  Niobium  Alloy 
(D-43) ,  Vanadium  Alloy(V-25  Cr-0.8  Zr) ,  and  Titanium  Alloy 
(Ti-6A1-4V) .  Properties  and  their  sources  are  summarized 
in  Appendix  A. 

2.  Limits  on  Data 

Properties  at  elevated  temperatures  are  available  on 
only  a  limited  basis.  Where  high  temperature  properties  are 
unavailable,  lower  temperature  values  are  used  as  the  best 
available  information.  Such  cases  are  indicated  by  the 
notes  in  Appendix  A.  Specific  high  temperature  data  will  be 
needed  and  may  be  incorporated  into  the  procedure  as  avail¬ 
able  to  achieve  a  more  precise  design  window  delineation. 

C.  Particle  Load  vs .  Lifetime  Tradeoff 

As  noted  in  Chapters  III  and  IV,  sputtering  and  fatigue 
pose  conflicting  constraints  on  the  allowable  material  thick¬ 
ness  for  any  geometric  configuration.  The  sputtering  erosion 
curves  of  Chapter  III  and  the  fatigue  limit  curves  of  Chapter 
IV  can  be  transferred  to  cycle  life  vs.  wall  thickness  axes 
for  various  heat  fluxes,  as  illustrated  by  figures  V-2  and 


V-3.  nc  that  fatigue  failure  is  unacceptable,  design 

conditions  must  lie  cn  or  below  both  the  sputtering  and 
fatigue  curves  for  a  specified  heat  flux.  The  intersection 
of  the  fatigue  and  sputtering  lines  is  seen  to  afford  the 
greatest  life.  Plotting  these  optimum  life  points  results 
in  the  design  curves  represented  by  the  solid  lines  of 
Figures  V-4  and  V-5.  The  behavior  of  material  optimum  life 
and  incident  heat  flux  is  apparent,  and  a  quick  determina¬ 
tion  of  the  target  capabilities  under  a  given  heat  flux  may 
be  made. 

In  Figures  V-2  and  V-3,  the  minimum  thickness  required 

to  contain  1000  psi  coolant  in  a  1  cm  inner  diameter  tube 

is  the  origin  of  all  sputtering  lines.  For  a  protective 

tile,  not  required  to  contain  pressure,  the  sputtering 

lines  would  pass  through  the  origin  and  have  the  same  slope 

as  in  Figures  V-2  and  V-3.  This  improves  the  optimum  life 

curves  slightly,  as  shown  by  the  dotted  curves  on  Figures 

V-4  and  V-5.  From  these  tile  curves,  the  optimum  material 

performance  limits  of  Figure  V-6  may  be  found.  Clearly,  the 

2 

goal  of  one  year's  survival  under  a  one  kw/cm  heat  flux  is 

not  met  by  any  of  the  materials.  Particle  delivered  heat 

2 

flux  would  have  to  be  reduced  to  .25-. 30  kw/cm  ,  (perhaps 
by  sloping  the  target  or  expanding  the  magnetic  field  lines \ 
before  a  one  year  survival  becomes  reasonable. 
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A _ k. 


D.  Single  Material  Analysis 

Simplicity  of  design  and  ease  of  manufacture  favor  a 

£ 

target  formed  of  a  single  material.  Composite  structures 
are  complicated  by  differential  thermal  expansion  and 
chemical  compatibility  considerations.  Accordingly,  analy¬ 
sis  of  a  simple  tube  and  a  simple  flat  plate  with  internal 
coolant  channels  is  first  undertaken. 

1.  Tube  Design 

Analysis  of  a  single  homogeneous  tube  is  achieved  by 
the  following  steps.  The  tube  is  considered  to  be  one  of 
many  in  a  uniform  array,  operating  under  the  base  case  eond 
itions  described  in  Chapter  1.  End  effects  are  initially 
neglected. 

a.  Thermo-hydraulic  design  limits  are  first  specified 
1)  The  operating  limit  due  to  critical  heat  flux 
is  established  using  the  modified  Lowdermilk  CHE  correla¬ 
tion  from  Chapter  II,  equation  7  .  (This  equation  uses 
British  Engineering  Units.) 


’crit 


r  _ . 05_  .15 

2 

.  „  T  (t+r ) 

D  L  q" 

4- 

4  *  L  hr ) 

90 

1400 

T  D  c 

t _ E _ _ 

(7) 


(T  -  T.  ) 
v  sat  in; 

2)  The  pumping  power  constraint  discussed  in 
section  II-C  is  applied. 


AfiL 


f  cp  *Ti-z 


£  .02 


till 


)  Clearest  visualization  of  the  hydraulic  design 
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window  is  afforded  from  a  plot  of  constraints  on  a  coolant 
temperature  rise,  vs.  channel  radius,  r,  set  of 

9 

axes.  From  equations  7  and  10,  the  critical  heat  flux 
and  pumping  power  ratio  limits  may  be  expressed  respective¬ 
ly,  in  British  Engineering  units,  as: 


CHF  : 


^crit-  '  V 


90  I" (S.F.)ct"  D*05L*15 
Gc  L  1400 


W  =.02  :  At, 


where ; 


Gc  l_  1400 

.637  k*357  L  q"*643 
p  r1.071  w  .357 
f  P 

.063  Ufb)*2  (/f)1,8 


fe) 


These  curves  are  plotted  for  base  case  conditions,  and 


(30a) 


appear  as  shown  in  Figure  V-7. 

b.  Thermo-mechanical  limits  for  each  material  are 
next  established. 

1)  The  minimum  wall  thickness  required  to  con¬ 
tain  internal  pressure  is  found  from  equation  20  as  a 
function  of  radius. 

2)  For  incremental  values  of  the  tube  radius,  r, 
the  thickness  corresponding  to  optimum  life  for  the  speci¬ 
fied  heat  load  is  derived  from  the  procedures  of  section 
V-C.  (The  goal  of  one  year  survival  must  here  be  reduced 
to  a  quest  for  longest  possible  life.)  Increasing  the  tube 
radius  raises  the  thickness  required  to  resist  hoop  stress. 


?h">rtons  the  o- 


life. 


t 


3)  An  upper  boundary  on  tube  radius  is  establish¬ 


ed  when  the  resulting  tube  life  falls  to  10  days. 

9 

4)  To  facilitate  manufacture  and  avoid  the  possi¬ 
bilities  of  bubble  blockage  in  the  tube,  a  minimum  tube 
radius  is  selected  as  .2  cm. 

5)  A  check  is  made  of  the  temperature  limits  at 
the  tube  exit  in  accordance  with  equations  27  and  28,  sec¬ 
tion  IV-E.  An  upper  bound  on  At  „  j.s  thereby  established, 

J.“4 

which  is  the  lesser  value  of: 


At 


1-2 


melt 


T  -T . 
melt  m 


m(l+t/r> 


q"/h 

(31) 


AT.  0  =  T  -T .  -  q"/h 

1-2 _ _ _  corr  m 

corr 


(32) 


c.  The  design  window  boundaries  are  plotted  for  each 
candidate  material,  as  illustrated  in  Figure  V-8a.  Materi¬ 
als  can  be  compared  by  plotting  the  window  for  fixed  life, 
as  in  Figure  V-8b.  This  constitutes,  the  simple  tube  design 
window.  Regions  within  the  window  are  acceptable  from  a 
hydraulic  and  mechanical  standpoint,  though  tube  life  de¬ 
creases  with  increasing  radius. 


2.  Plate  Design 

Assuming  that  plates  are  restrained  from  bowing  but 
can  expand  in  directions  parallel  to  their  surface,  essen¬ 
tially  the  same  analysis  can  be  performed  for  the  flat 
plate  as  was  performed  for  the  tube.  An  assumption  of 
rectangular  coolant  channels  would  require  greater  wall 


thicknesses  to  resist  internal  pressure.  Circular  channels 
may  result  in  added  thermal  stress  in  the  thick  mat^  rial  be¬ 
tween  tribes.  This  increased  stress  could  be  relieved  by 
sloped  grooves  in  the  face  of  the  plate,  so  that  the  tube 
analysis  remains  an  appropriate  procedure.  In  this  config¬ 
uration,  the  spacing  between  channels  is  independent  of  the 
thickness  facing  the  plasma.  The  channel  spacing  may  be 
selected  to  optimize  the  conduction  to  the  rear  of  the 
channel,  reducing  pumping  power  requirements  per  unit  area 
of  target.  Since  sputtering  and  fatigue  will  affect  the  plas 
ma  facing  wall  in  the  same  way  as  the  tube,  the  mechanical 
temperature  and  radius  limits  will  be  identical  to  the  tube 
case.  The  design  window  would  appear  as  in  Figures V-8. 

E.  Composite  Material  Analysis 

The  diverse  demands  for  material  resistance  to  sputter¬ 
ing,  fatigue,  and  erosion,  and  for  compatibility  with  plasma 
and  coolant,  suggest  that  some  blend  of  materials  may  best 
suit  the  divertor  purpose.  As  seen  in  the  previous  section, 
no  single  material  offers  adequate  properties  for  all  con¬ 
straints  in  the  base  case  design. 

Analysis  of  composite  designs  encompasses  all  the  con¬ 
siderations  of  the  single  material  targets.  In  addition,  the 
effects  at  interfaces  between  different  materials  must  be 
assessed. 


1.  Coated  Tube  Analysis 

Based  on  the  results  of  section  V-D,  the  concept  of 


adjacent  -<-as  seen  to  be  somewhat  less  desirable  than 

circular  channels  in  a  flat  plate,  since  the  thickness  of 
tube  walls  to  resist  sputtering  was  linked  to  the  spacing 
between  coolant  channels.  Very  thick  tubes  would  result  in 
an  increased  heat  load  per  channel.  However,  if  sputtering 
resistance  could  be  provided  by  an  appropriate  coating  on  the 
plasma  side  of  a  tube  array,  the  heat  load  per  tube  could  be 
fixed,  and  a  coated  tube  array  would  be  competitive  with  flat 
plate  concepts  by  virtue  of  its  production  simplicity. 

Thermohydraulic  analysis  would  be  identical  to  that  of 
the  simple  tube  case.  The  minimum  wall  thickness  for  the 
tube  substrate  is  determined  as  before,  based  on  the  hoop 
stress  consideration.  The  coating  thickness  must  be  given 
particular  attention,  in  order  that  it  neither  crumbles  due 
to  fatigue  nor  is  penetrated  by  sputtering.  Norem  and  Bowers 

(47)  report  that  a  thickness  of  just  10  microns  of  beryllium 
should  be  sufficient  to  accomodate  incident  ions  without 
sputtering  from  the  substrate.  They  report  that  such  a  thin 
coating  would  be  redeposited  every  20  pulses,  assuming  ade¬ 
quate  coating  material  is  artificially  introduced  into  the 
plasma  edge.  Such  a  thin  coating  should  not  be  subject  to 
thermal  stress,  since  its  surface  would  be  rough  on  a  scale 
of  Angstroms  (Ref.  47).  The  development  of  this  sort  of  in- 
situ  plasma-recoating  of  the  target  surface  would  virtually 
eliminate  the  sputtering  limits,  making  the  coated  tube 
more  attracimv^. 
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2.  Armored  Plate  Design 

Thermohy  dr aul  i  c  analysis  of  the  armored  plate  would  be 

9 

identical  to  that  of  the  simple  flat  plate.  Thermo-mechanic- 
ally,  procedures  would  parallel  the  simple  case,  except  that 
surface  stresses  arising  at  the  interface  between  the  two 
materials  must  be  accounted  for.  If  the  armored  tile  is 
bonded  to  the  substrate,  the  analysis  illustrated  by  Figure 
V-9  affords  a  two  dimensional  indication  of  the  resulting 
thermal  stresses.  For  the  substrate,  the  temperature  differ¬ 
ence  through  the  thickness  is: 

AT1  =  ^k-*1  *  <33) 

Likewise,  for  the  tile: 

q"  t2 

fiT,  =  - — .  (34) 

^  k 

*2 

Were  the  plates  unrestrained,  they  would  appear  as  in  Figure 
V-9c  under  a  steady  heat  load.  Imposing  a  plane  strain 
condition  on  the  plate  leads  to  a  virtual  thermal  strain  as 
shown.  The  neutral  plane  through  the  plate  will  have  dis¬ 
placed  a  distance  X  from  its  unheated  position  on  the  x  axis 
due  to  thermal  expansion.  The  virtual  strain  in  the  substrate 
is  then: 

<1  =  X  -  ATj.  ^  ATX  )  ,  (35) 

while  that  in  the  tile  is: 

<2  -  X  -  («2  4Tj  .  (361 

For  the  two  dimensional  plate,  the  stress  in  each  region  is 

t:  V  43): 
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,  C'1  =  1 


T  2  = 


- i [x  -  (*  AT  ~  +  «  AT  )1  (37) 

1  -  Vj )  A  1  ^1  1  A  J 

K 

- ? Tx  -  (*  AT_-2-  +  «  AT.)1  (38) 

( 1  -  V,  )  L  2  z  1  J 


For  equilibrium  we  can  set  the  sum  of  forces  on  the  plate 
equal  to  zero  and  solve  for  the  displacement  X: 

’o  E, 


Force  =  0  = 


(39) 


-t,  ^l 


'2  E, 


1“ 


E,  t- £.T, 

(1-1),)  i-  M— ^ 


(40)  X  = 


"1AT1  1 

-  VT1  -  ^  2J  < 

.-"-t-J 


X  -  o<2At 


dz 


’E_t,c<  I 

"2-- - -  +  E2t2^2ATj  (1-tf.) 


(1  -  ))2)  E1t1  +  E2t2(l  - 


This  value  of  X  may  then  be  used  to  determine  the  strains  and 
stresses  from  equations  35,  36,  37,  and  38.  Sample  results 
appear  as  in  Figure  V-10.  The  peak  strains  so  determined  can 
be  compared  with  the  maximum  allowable  strain  range  for  the 
specified  material  and  design  life.  Conversely,  the  maximum 
allowable  strain  may  be  used  in  equations  33,  34,  35,  and 
36  ,  to  determine  the  maximum  allowable  heat  flux  for  a  spec¬ 
ified  target  thickness.  In  this  way,  fatigue  curves  on  a 
life  vs.  thickness  axis  may  be  generated,  the  sputtering  rate 
plotted,  and  optimum  life  vs.  thickness  curves  obtained. 


Analysis  of  six  types  of  tiles  (Aluminum,  TZM,  Beryllium, 
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Graphite,  Titanium,  and  Vanadium)  over  a  copper  substrate 
were  analyzed  in  this  way.  In  all  cases  except  aluminum,  the 
tensile  virtual  strains  in  the  tile  were  reduced  over  the 
unbonded  plate.  The  compressive  strains  increased,  but  oc¬ 
curred  in  the  first  layers  of  surface  to  be  eroded  by  sput¬ 
tering.  A  modest  improvement  in  the  tile  fatigue  life  should 
result  from  this  strain  reduction.  The  relatively  high 
expansion  coefficient  of  copper  tends  to  mitigate  the  back 
side  tensile  strain  of  all  candidate  tiles  except  aluminum, 
reducing  cyclic  fatigue  at  the  back  of  the  tile.  With  remote 
replacement  of  tiles,  the  substrate  and  coolant  channels 
would  easily  survive  past  the  one  year  design  goal.  The 
relative  ranking  of  material  survival  would  be  unchanged 
from  Figure  V-6,  though  life  would  be  increased  slightly. 

3.  Mechanically  Unbonded  Layered  Plate  Design 

The  foregoing  analysis  has  revealed  serious  limitations 
on  the  lifetime  and  heat  load  constraints  under  which  solid 
target  configurations  may  operate.  An  upper  limit  on  the 
ability  of  a  solid  target  to  survive  may  be  evaluated  by 
assuming  a  design  which  in  effect  takes  full  advantage  of 
material  temperature,  sputter  rate,  and  fatigue  thickness 
constraints  simultaneously.  This  is  achieved  by  assembling 
relatively  thin  layers  of  armor  material  so  that  mechanical 
bonding  between  them  is  non-existent,  while  thermal  bonding 
remains  very  good.  Here  a  thin  film  of  liquid  lithium  is 
pc. .relate  !  o  separate  the  layers.  Powdered  graphite  may  be 


m  accoy .  •  i  torn  uti  vu.  layer  is  thin  enough  to  pre¬ 

vent  fatigue  failure  during  its  survival  time  in  the  target. 

.  9 

In  this  way  successively  thicker  layers  can  be  built  upon 
the  cooled  backing,  up  to  a  thickness  at  which  a  temperature 
limit  is  approached.  For  most  materials  considered,  this 
limit  is  the  boiling  point  of  lithium.  The  thickness  of 
armor  material  thus  assembled  resists  sputtering  for  the 
longest  time  possible  for  any  target  configuration. 

The  obvious  price  paid  to  reach  this  optimum  life 
design  is  an  increase  in  the  complexity  and  difficulty  in 
fabrication  of  the  target.  To  effectively  eliminate 
mechanical  bonding  between  solid  layers,  the  liquid  lithium 
would  require  an  expansion  reservoir  connected  to  each  liquid 
layer.  These  connections  may  have  to  be  valved  in  some  way 
to  prevent  loss  of  lithium  as  sputtering  penetrates  to  the 
depth  of  each  lithium  layer.  Thus  as  few  layers  of  lithium 
as  possible  would  be  desirable  to  simplify  valve  assembly. 

Figure  V-ll  indicates  the  analysis  of  this  optimum 

configuration  for  the  eight  candidate  materials  of  this 

2 

report.  TZM  affords  the  longest  life  at  a  1  kw/cm  heat 
load;  261  days.  Niobium  (D-43)  appears  to  be  nearly  as 
attractive  at  258  days,  perhaps  more  so  since  but  four  lith¬ 
ium  layers  are  required  as  opposed  to  TZM's  nine  layers. 

If  an  ideal  lubricant  between  the  layers  (k  =  T 
=  «*>)  were  available  the  total  wall  thickness  could  be  in- 

..ur  cner ,  up  u  uu  ..men  the  oacei  auieucs 


This  would  improve 


approaches  the  material  melting  point, 
the  target  lifetimes  by  roughly  a  factor  of  2  over  the  values 
of  Figure  V-ll.  This  extreme  would  probably  be  impractical, 
however,  since  the  remainder  of  the  divertor  duct  would 
become  even  hotter  due  to  radiative  heat  transfer. 
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SENSITIVITY  ANALYSIS 


$ 

The  base  case  parameters  selected  for  the  initial 
analysis  in  Chapter  I  are  typical  of  the  conditions  current¬ 
ly  being  considered  for  divertor  operation.  The  resulting 
design  window  provides  a  view  of  how  much  option  space 
exists  for  these  conditions.  Much  latitude  exists  for 
varying  many  of  these  parameters,  however,  in  an  effort  to 
expand  the  design  window  and  determine  its  behavior  over 
a  broad  range  of  conditions.  Further,  uncertainty  in  the 
applicability  of  presently  known  physical  behavior  of  cool¬ 
ants  and  structural  solids  at  the  extreme  conditions 
anticipated  suggests  that  better  knowledge  of  such  behavior 
may  affect  the  design  window  boundaries.  Accordingly,  a 
systematic  analysis  of  the  effect  of  varying  key  parameters 
and  physical  assumptions  has  been  made. 

A.  Pulse  Length  and  Fatigue  Constraint 

The  most  severe  limitation  on  solid  divertor  target 
options  at  high  heat  loads  appears  to  be  the  dual  constraint 
of  sputtering  and  cyclic  fatigue.  Extension  of  the  reactor 
pulse  length  would  be  a  means  of  easing  the  fatigue  problem. 
An  examination  of  the  improvement  in  material  lifetimes 
expected  from  an  order  of  magnitude  reduction  in  the 
number  of  thermal  stress  cycles  imposed  on  the  target  per 
ye^r  was  performed.  Allowable  strain  ranges  for  fatigue 
lives  between  10J  and  104  cycles  from  Ref.  42,  combined 
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v.'j  f:  the  equation, 
appropriate  fatigue 


t  < 


2  k 


life  vs. 


- - ^ ,  established  the 

ex  q 

thickness  curves  for  various 


heat  fluxes. 

Sputtering  rates  at  each  heat  flux  were  determined  as 
in  section  III-A.  At  each  heat  flux  the  thickness  and 
corresponding  lifetime  at  which  the  minimum  thickness  to 
withstand  sputtering  just  equaled  the  maximum  thickness  to 
withstand  fatigue  was  the  optimum  material  design.  For 
niobium,  these  points  appear  or.  Figure  VI-1  along  with  the 
base  case  optimum  life  curve.  For  a  given  heat  flux,  an 
order  of  magnitude  reduction  in  cycles  per  year  provides 
only  80%  increase  in  the  target  lifetime.  Improvements  in 
other  materials  were  of  the  same  order. 


B .  Sputtering  Characteristics 

The  devastating  rate  at  which  target  material  appears 
to  be  eroded  by  sputtering  in  this  analysis  demands  a 
search  for  ways  in  which  the  situation  may  be  improved. 
Several  factors  may  be  considered.  First,  there  is  some 
scatter  in  the  reported  data  for  sputtering  coefficients. 

A  comparison  of  data  from  Ref.  33  with  corresponding  predic¬ 
tions  of  the  analytical  formula  of  Ref.  51,  Figure  VI-2, 
indicates  a  factor  of  two  uncertainty  in  these  coefficents. 
Second,  the  imposed  plasma  load  may  be  made  less  severe  than 
is  currently  predicted  by  the  INTOR  study.  Since  light  ion 
sputtering  coefficients  peak  near  energies  of  1  keV, 

.  ..-..e  .i  j  v.'  icing  the  ion  energies  woulc.  reduce  ciieir 


sputter  ir.:.  -c  ••  -'f  icient,  thereby  reducing  the  erosion 
rate.  At  the  same  incident  heat  flux,  fewer  particles  of 
higher  ehergy  or  more  particles  of  lower  energy  could  re¬ 
duce  the  sputtering  rate  by  up  to  50%.  Taking  as  an 
optimistic  improvement  a  four- fold  reduction  in  the 
sputtering  rate,  the  more  favorable  optimum  life  curve  of 
Figure  VI- 3  results.  Lifetimes  for  each  heat  flux  are  seen 
to  increase  approximately  2.5  times. 

Such  optimism  may  be  warranted  in  light  of  suggestions 
that  radiation  in  the  plasma  might  be  enhanced,  resulting  in 
less  energetic  particles  at  the  divertor  target.  (Ref. 13) 
However,  the  degree  to  which  relatively  high  Z  materials  may 
be  ionized  and  accelerated  into  the  divertor  target  will 
adversely  affect  the  sputter  rate.  Indeed,  ions  of  Z >  6  can 
achieve  sputtering  coefficients  greater  than  one  if  accele¬ 
rated  to  1  keV,  possibly  causing  a  catastrophic  sputtering 
cascade.  Finally,  the  divertor  target  may  occasionally 
receive  the  intense  particle  load  of  a  plasma  disruption. 

The  conditions  in  the  divertor  chamber  affecting  self-sput¬ 
tering  and  plasma  disruption  effects  are  ill-defined,  and 
have  not  been  included  in  this  analysis. 

C.  CHF  Limit 

Due  to  present  uncertainties  in  the  prediction  of  sub¬ 
cooled  boiling  pressure  drop,  the  critical  heat  flux  limit 
for  the  design  window  has  been  restricted  such  that  the 
pieuxcted  critical  heat  flux  is  1.3  times  larger  tnan  the 
design  heat  flux.  (See  section  II-B.)  Should  experience 


with  subcooled  liquids  in  high  heat  fluxes  provide  confidence 
that  this  safety  factor  nay  be  reduced,  an  appreciable  enlarge¬ 
ment  of  the  design  window  operating  region  will  result.  Figure 
VI-4  indicates  the  effect  of  varying  the  safety  factor  from  1.1 
to  1.5.  A  wider  selection  of  operating  conditions  becomes  a- 
vailable  as  safety  factor  is  lowered.  Specifically,  reducing 
the  safety  factor  from  1.3  to  1.1  while  maintaining  other  base 
conditions  reduces  the  required  mass  flux  by  30%. 

The  prospect  of  operating  a  high  mass  flux  system  under 
conditions  of  subcooled  boiling  on  the  tube  walls  introduces 
the  potentially  catastrophic  occurrence  of  pressure  drop-flow 
rate  instabilities.  Should  the  flow  rate  decrease  momentarily 
due  to  any  perturbation,  the  increased  boiling  would  increase 
pressure  drop,  in  turn  lowering  Tgat  and  possibly  driving  the 
tube  to  a  burnout  condition.  Lowdermilk,  et  al,  among  others 
dealt  extensively  with  this  phenomenon  in  their  small  tube, 
high  heat  flux  experiments.  (Ref.  20.'  Throttling  the  inlet 
to  each  heated  tube  would  help  control  such  instabilities,  but 
would  incur  a  significant  pumping  power  penalty. 


D.  Pumping  Power  Limit 

Relaxation  of  the  pumping  power  ratio  limit  from  2%  to 

5%  or  10%  may  prove  necessary  and  justifiable,  particularly  in 

2 

cases  where  heat  fluxes  greater  than  1  kw/cm  may  be  received. 
The  analysis  here  does  not  account  for  pressure  losses  in 
headers,  bends,  throttles,  etc.,  so  pushing  the  pumping  power 
limit  ajjove  20%  would  seriously  impact  on  energy  economy. 


Such  a  limit  relaxation  would  clearly  widen  the  operating 
region  of  the  design  window,  allowing  a  greater  margin  of 
safety  against  CHF,  as  shown  by  Curves  A  and  B  in  Figure  VI-9. 

E.  Inlet  Temperature 

Increasing  the  inlet  temperature  decreases  the  exit 

subcooling.  This  necessitates  higher  mass  flux  to  avoid  CHF. 

The  accompanying  increase  in  pumping  power  pushes  the  CHF 

constraint  ever  closer  to  the  W  =  .02  constraint.  With 

P 

other  parameters  unchanged  from  the  base  case,  the  design 
window  vanishes  for  T^n  larger  than  190°C  Operation  at  an 
inlet  temperature  as  close  to  freezing  as  possible  widens  the 
design  window  and  allows  operation  at  the  lowest  possible 
mass  flux,  as  evidenced  by  Figure  VI-5. 

F.  Pressure 

Varying  pressure  changes  the  degree  of  subcooling  for  a 
fixed  inlet  temperature..  High  pressures  provide  an  increase  in 
subcooling,  allowing  lower  mass  flux  to  avoid  CHF.  This  en¬ 
larges  the  hydraulic  design  window,  as  shown  in  Figure  VI-6, 
but  adds  to  the  hoop  stress  imposed  on  the  material.  The 
radius  limit  due  to  hoop  stress  imposed  on  the  structure  is 
proportionately  reduced.  Further,  the  occurrence  of  system 
pressure  cycling  would  impose  an  additional  life-shortening 
fatigue  component  to  the  structural  analysis.  Finally,  a 
preference  for  low  to  moderate  pressure  systems  is  indicated 
by  the  consequenc »3  of  a  high  pressure  pipe  break  within  the 
divertor  chamber. 


G.  Tube  Length 

Variation  of  the  tube  length,  keeping  other  parameters 

« 

constant,  was  investigated.  The  resulting  operating  regions 
are  illustrated  by  Figure  VI-7  for  lengths  of  10cm,  50cm,  and 
100cm.  Because  the  temperature  rise  in  the  tube  is  greater 
for  longer  lengths,  less  exit  subcooling  remains.  A  higher 
mass  flux  is  necessary  to  avoid  CHF.  Entrance  effects  become 
less  significant  at  longer  lengths,  improving  the  pumping 
power  ratio  over  the  base  case.  The  net  effect  is  that  the 
permissible  operating  region  shrinks  as  length  is  increased. 

H.  Thickness 

Due  to  the  wall  thickness  demands  of  sputtering  and 
fatigue,  adjacent  tube  arrays  will  be  required  to  operate  with 
considerable  wall  thickness  separating  the  coolant  channels. 

As  this  thickness  increases,  the  linear  power  received  by  each 
channel  also  increases;  9' tot  =  ^"act  *  However,  the 

conduction  path  to  the  rear  side  of  the  channel  is  improved, 
permitting  extraction  of  more  heat  by  the  coolant  at  essentially 
the  same  peak  heat  flux  on  the  hottest  side.  The  pumping  power 
ratio  is  thus  improved,  resulting  in  an  enlargement  of  the  de-" 
sign  window  for  wall  thicknesses  between  ,05(r)  and  3.0(r). 

For  wall  thicknesses  above  3.0(r),  heat  flux  to  the  coolant 
exceeds  the  incident  heat  flux  from  the  plasma,  q"  ,  raising 
mass  flux  and  shrinking  the  design  window.  An  optimum  thick- 
nr.~fr  -1* y  be  determined  car  ach  structural  material  at  which 
the  design  window  is  at  its  "widest".  See  Figure  VI-8. 


I.  Heat  Flux 


An  investigation  of  the  limits  to  which  steady  state 
heat  flux  may  be  raised  was  undertaken,  keeping  all  other 
parameters  as  in  the  base  case.  Extrapolating  the  CHF  cor¬ 
relation  based  on  Lowdermilk's  work  to  higher  heat  fluxes 
imposes  a  very  stringent  constraint  on  the  design  window. 

As  shown  in  Figure  VI- 9,  the  increased  mass  flux  required  to 
accomodate  a  4  kw/cm  heat  flux  raises  the  minimum  pumping 
power  ratio  to  avoid  CHF  to  approximately  60%.  At  these 
higher  heat  fluxes,  the  correlation  on  Rousar  (Ref.  19)  is 
more  appropriate  for  CHF,  and  leads  to  the  operating  regions 
depicted  in  Figure  VI-10.  Clearly,  a  significant  pumping 
power  penalty  is  incurred  by  raising  heat  flux.  The  work  of 

M.A.  Hoffman,  et  al.,  (Ref.  12),  provides  an  estimate  of  the 

.  2 
option  space  available  at  a  3.3  kw/cm  heat  flux,  and  is  in 

close  agreement  with  Figure  VI-10. 
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Figure  VI-1:  Sensitivity  to  Change  in  Pulse  Duration 
of  Optimum  Life  vs.  Thickness  Curves 

;or  Niobium  Target. 
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Figure  VI-2:  Comparison  of  Experimental  Data  with 
Analytical  Results  for  Sputtering 
Coefficient. 
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VII.  CONCLUSIONS 


9 

The  conclusions  of  this  investigation  of  solid  divertor 
target  options  can  be  summarized  as  follows. 

1.  A  hydraulic  system  using  non-boiling  highly  sub¬ 
cooled  water  appears  feasible  to  cool  surface  heat  loads  of 

2 

1  kw/cm  in  the  divertor  target  application.  Intermediate 
pressures  (200  to  500  psia)  are  required  to  provide  adequate 
subcooling  to  avoid  nucleate  boiling.  Channel  lengths  from 
10  cm  to  100  cm  may  be  used. 

2.  Provision  for  handling  hot  spots  and  off-normal 

2 

conditions  with  heat  loads  up  to  2  kw/cm  may  be  made  by 
operating  the  system  with  very  high  mass  flux  and  high  inlet 
subcooling.  The  pumping  power  to  heat  transfer  ratio  can 
be  kept  at  2%.  Careful  design  will  be  necessary  to  avoid 
internal  cavitation  erosion. 

3.  The  Saha-Zuber  equation  for  "net  vapor  generation" 
compares  favorably  with  existing  CHF  correlations  in  the 
range  of  variables  pertinent  to  divertor  applications.  It 
provides  a  theoretical  basis  for  predicting  CHF  and  should 
prove  valuable  in  extending  CHF  prediction  beyond  the  ranges 
of  existing  data. 

4.  At  high  mass  flux,  the  onset  of  fully  developed 
nucleate  boiling  will  increase  the  total  pressure  drop  over 
the  non-boiling  case  two  to  six  times.  This  can  be  accomo- 

uauu  in  channels  oj.  u-La.tieter  greater  chan  4  nun  uiid  ac 


pressures  of  500  psia  with  less  than  a  10%  overall  pressure- 
drop.  The  likelihood  of  ?  pressure  drop-flow  rate  instabil¬ 
ity  in?  a  bank  of  flow  channels,  however,  remains  an  obstacle 
to  operation  in  this  regime. 

5.  The  cyclic  thermal  loading  and  ion  sputtering 
conditions  currently  envisioned  at  the  divertor  target 
surface  limit  the  life  of  conventional  tube  or  plate  targets 
for  the  eight  candidate  materials  studied,  to  several  months 
Molybdenum  Alloy  (T2M) ,  and  Niobium  Alloy  (D-43)  appear  most 
promising  from  a  materials  standpoint,  but  may  well  quench 
the  plasma  unless  ions  sputtered  from  the  target  are  effect¬ 
ively  prevented  from  penetrating  the  plasma  core.  Beryllium 
appears  most  favorable  among  the  low  Z  candidates,  but  can 
survive  only  one  month  without  redeposition. 

6.  Development  of  a  mechanically  unbonded,  thermally 
bonded  laminated  design  would  allow  target  life  to  be  ex¬ 
tended  to  9  months. 

7.  A  copper  substrate  protected  from  sputtering  could 
be  designed  to  survive  the  cyclic  loading  for  well  over  one 
year. 

8.  Techniques  for  rapid  remote  replacement  of  the 
divertor  target,  in-situ  recoating  of  the  surface,  reduction 
of  the  particle  load,  or  steady  state  tokamak  plasma  opera-' 
tion  appear  essential  to  development  of  a  solid  divertor 
target  capable  of  surviving  for  one  year  without  plant 


VIII.  RECOMMENDATIONS 


The  investigation  of  divertor  target  design  and  cooling 
herein  described  has  revealed  numerous  areas  in  which  fur¬ 
ther  effort  is  necessary.  The  limits  on  lifetime  and  heat 
load  predicted  by  this  general  approach  impose  severe  re¬ 
strictions  on  plant  capacity  factor,  since  with  conventional 
tube  or  flat  plate  arrangements  frequent  plant  shutdowns  may 
be  necessary  to  allow  refurbishment.  Methods  of  in-situ 
coating  of  the  divertor  target  during  rejuvenation  periods, 
if  developed,  might  greatly  ease  design  by  alleviating  the 
wall  thickness  constraint  due  to  sputtering.  Without  such  a 
technique,  the  development  of  a  laminated,  thermally  bonded, 
mechanically  unbonded,  actively  cooled  flat  plate  target 
should  be  emphasized  as  a  promising  solid  target  option. 
Targets  mounted  on  central  supports,  with  edges  free  to 
bend  and  expand,  should  be  examined  as  a  possible  means  of 
reducing  thermal  fatigue. 

Areas  of  investigation  necessary  for  improved  evalua¬ 
tion  of  all  solid  target  options  include  the  following; 

1.  Validity  of  Saha-Zuber  prediction  of  "net  vapor 
generation"  as  an  accurate  forecast  of  CHF  conditions  over 
a  broad  range  of  variables. 

2.  Effect  of  non-uniform  temperature  distribution  in 
coolant  on  fluid  viscosity  and  overall  pressure  drop. 

3.  L.ic.er standing  and  correlation  of  subcooled  boiling 


pressure  drop  in  coolant  channels. 

4.  Thermal  stresses  at  tube  connections  and  plate 
edges . 9 

5.  Material  properties  at  elevated  temperatures, 
particularly  modulus  of  elasticity. 

6.  Low-cycle  fatigue  life  data  verification  for 
candidate  materials  at  elevated  temperatures. 

7.  Occurrence  and  behavior  of  redeposition  after 
sputtering. 

8.  Aggregate  effects  of  simultaneous  sputtering  and 
cyclic  fatigue  under  surface  heat  loads.  (Does  sputtering 
promote  or  retard  fatigue  damage?) 

9.  Chemical  erosion  of  graphite  exposed  to  hydrogen 
plasma . 

10.  Transport  of  sputtered  particles  from  target  to 
plasma  core. 

11.  Distribution  of  particle  fluxes  and  energies  across 
the  target  surface. 

12.  Techniques  for  helium,  hydrogen,  and  impurity 
separation. 

Areas  of  investigation  specific  to  the  suggested 
laminated  target  design  are  listed  below: 

13.  Effect  of  released  lithium  on  plasma  burn. 

14.  Control  of  tritium  generated  in  lithium  layers. 

15.  Positive  closure  valves  for  lithium  expansion  ducts 

16.  Mechanical  behavior  at  constraining  edges. 
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PPENDIX  A 


::h  cals  properties 


Table  ,-\-l  provides  representative  property  values  for 
6 

the  eight  candidate  materials  in  the  range  of  operating 
temperatures  indicated.  Physical  properties  are  evaluated 
at  the  midpoint  of  the  temperature  range  for  each  material, 
except  that  strain  and  stress  limits  are  taken  so  as  to 
give  a  conservative  bound  over  the  entire  operating 
temperature  range.  The  improvement  in  material  lifetime  a- 
va liable  from  abandoning  this  conservatism  is  less  than  a 


factor  of  two. 

friaterials  data  is  taxen  from  References  3o*  42,  43, 


44,  49,  and  50.  Pertinent  data  for  water  and  lithium  is 
taken  from  References  23  and  57* 

The  development  of  the  fatigue  life  versus  thickness 


curves  appearing  in  Figures  V-2  and  V-3  utilises  equation 


26  and  the  maximum  strain  range  data  for  fatigue  life. 


t<  z  *  O-h 


q* 


max 


(26) 


As  an  example,  substituting  in  property  values  for  Nb  and 

q"=l  kw/cm2  yields  t  5  [lof  ^  max  =  * 

The  following  data  points  result. 


Ufe 

10**  cycles 
2x10^  cycles 
4x10^  cycles 
10*  cycles 


*C.(0 

.210 

.180 

.160 


1*9 

1.5 

1.4 


Plotting  these  points  on  the  life  vs.  thickness  axes  yields 

o 

the  fatigue  curve  for  1  kw/cm  of  Figure  V-3.  Changing  q" 


Jl  i  t  Vl  Ofi  OVlJv1  J 


w  u  «.  *  1  or  C*  ~  -  C*- 


fluxes 


Operatic?.  Temp.- 
Range  ”(°C) 

Water  Corrosion 
Temp.  Limit  (°C) 
Mieit  Temperature 

k  (w/m-°K) 

E  (CPa ) 

(xl0~°/  :;’K) 

Poisson* 3  Ratio,  9 

Density,  p 
(g/cm3) 

Atomic  Weight,  M 

Allowable  Stress 
Intensity  (MFa) 

Max  Strain  Range  (%) 
for  Fatigue  Life: 
10^  cycles 
2x10^  cycles 
4x10^  cycles 
5x10^  cycles 
10->  cycles 


200- 

800 

300 

1290 


400- 

2000 


200- 

400 


300- 

400 

360 


200- 

800 

360 


300-| 

600 


200- 

800 

330 


200- 

1000 

360 


1668  1900  1085  2468  2610 


15.0  10.0  25.0  10.2  10.2  18.0 


1.85  2.10  2.70  4.54  6.11  8.96  8.57  10.2 
9.01  12.0  27.0  47.9  50.9  63.5  92.9  95.9 


.116 


|  .^91 

1  .  * 


.083  *063  .240  .614  .343  .278  .070 

.073  .226  .290  .210 

.080  .220  .240  .180 

.421 

.061  .089  .215  .376  .200  .160  .046 

♦  * 


*  indicates  room  temperature  fatigue  data;  high 
temperature  data  not  available. 

Sputtering  coefficients  provided  in  Table  III-l. 


Table  A-l 1  Representative  Materials  Properties 


Heat  Flux,  qM 

10  Mw/m2 

1  kw/cm^ 

3.17x1 06  BTU 
hr-f t^ 
7.37xl06  lbm 

Mass  Flux,  G 

io4  kg 

10.3  g 

sec-m * 

sec-cm^ 

hr-ft^ 

Mass  Flow  Hate,  m 

1  kg/sec 

103  kg/sec 

7920  lbm/hr 

Volume  Flow  Rate 

472  cm3 
sec 

7.48  GPM 

1  ft3/min 

Pressure,  p 

.1  MPa 

\?50  torr 

14.5  psia 

Length,  L 

1  m 

100  cm 

3.281  ft 

tess ,  m 

1  leg 

1000  g 

2.2046  lbm 

Power,  P 

1  Mw 

1000  kw 

3. 41x1 06  BTU 
hr 

Density, p 

3  g/cm3 

1000  kg/m3 

62.4  lbm/f t3 

Therrral  Conductivity, 

100  w 

.239  cal 

57.8  .  BTU 

k 

ro^C 

sec-cm-u< 

D  hr-ft-uF 

Specific  Heat,  cp 

1  eal/g-°C 

s 

| 

1  BTU/lbm-°F 

Viscosity, M 

10  poise 

|l  kg/sec -m 

2419  lbm/hr-1 

Table  A-2«  Conversion  Equivalents 
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APPENDIX  Bt  fcETHODOJX-ay  oUi  ’vURY 


^nis  appendix  summarizes  the  design  window  methodology 
used  in  this  investigation.  The  primary  analysis  is  done 
for  parallel  tube  arrays  exposed  to  particle  ana  heat  loads 
on  one  side  only.  The  procedure  is  then  modified  to  treat 
protective  tiles  and  flat  plates  with  internal  channels. 

The  thermo-hydraulic  constraints  are  first  determined. 
Tube  lengtn  and  system  pressure  are  selected  and  heat  load 
is  specifier,'.. 

1.  Critical  heat  flux  is  determined  using  the  Lowder- 
milk  correlation  (Ref.  20)  modified  to  account  for  varied 
subcooling,  after  Rousar  (Ref.  19) •  A  safety  factor  of  1.3 


is  applied  to  CKF. 


fcrit 


d,q5l,: 

/O' 


^  £ 

q"  '2  q"  L 

+  D 

(^sat  "  ^in^ 


(7a) 


*  Note  1  The  procedure  of  Saha  and  Zuber  may  be  substituted 
here  to  predict  CHF  on  a  more  theoretical  basis.  See  Ref. 54. 
2.  Functional  dependencies  among  variables  are  estab¬ 


lished  by  conservation  of  energy  in  the  tube. 


ATI-* 


4  q"  L  (t+r)/r 
IT  Q  D  cn 
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3.  The  pumping  power  constraint  is  next  established. 

-  *  .02  (11) 


WD  «  - £E— -  *  .02  (11) 

P  J>  cp  atx.2 

The  limit  of  2 %  is  a  soft  constraint,  and  may  be  exceeded 


/  fo;  successful  divertor  operation. 
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h.  The  existence  of  subcode^  coiling  in  the  tube  is 

eva lusted  using  the  Rohsenow-Bergles  correlation  for  onset 

of  nucleate  boiling  in  water. 
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1  15.nO  pl«i56  J 

5.  If  G  is  greater  than  the  right  hand  side  of  equa¬ 
tion  1,  the  pressure  drop,  Ap,  is  determined  using  the  con¬ 
ventional  barcy-A/eisbach  approach  with  a  viscosity  correc¬ 
tion  for  non-isothermal  flow. 

APnonb  -  fj-  Cl.  ■  »  Sl± 

2  D  p  gc  2  D1,2  j>  gc 

*  Note*  The  nature  of  this  viscosity  correction  factor 
needs  further  refinement  for  the  non-uniform  temperature 
conditions  found  in  the  divertor  target.  See  Ref.  56. 

6.  If  G  is  less  than  the  right  hand  side  of  equation 
1,  the  location  of  ONB  must  be  determined  and  the  Mendler 
procedure  for  determining  subcooled  boiling  £p  employed 
over  the  subcooled  boiling  length. 


^  Pscb  “ 


2  6, 


D  J Lj  d" 


(8a) 


with! 


1—  *  (  1-.0025  ATj_l)  Ql  +  ,76(G/106) 


-2/3 


iso 


v  1  — - 


AT 


J-L 


n 


•  Voovq "/n)  'J 


(  3b  j 


1  ?'} 


ATj-l 


sax 


■  m 


+  rc  (qyio6)^  {tic) 
vr  G  Cp  D  exp( p/900) 


is  given  in  equation  2. 

?.  Equations  ?j  and  11  constitute  the  principal 
thermohydraulic  constraints. 

8.  Thermo-mechanioal  constraints  are  next  imposed. 
The  sputtering  erosion  rate  is  determined. 

-A-t_  =  t  SiJjC  (-L- )  U3) 

<  Vna 

9*  Adding  a  hoop  stress  constraint  gives  a  minimum 
wall  thickness  requirement. 

*  -  sSt-  +iTsiJic  (A>  (l3a) 


10.  The  fatigue  constraint  is  found  from  the  following 
equation,  where  AC max  varies  with  target  life. 

t<  -i-Ulli-L  A£.ax  (26) 

tx  q" 

11.  Solving  13a  and  26  simultaneously  for  the  design 
heat  load  determines  the  optimum  life  and  wall  thickness  for 
the  material.  Varying  radius  changes  the  maximum  life  of 
the  material.  Fixing  life  provides  an  upper  limit  on  tube 
radius. 

12.  Temperature  limits  are  checked.  These  may  pose 
an  upper  limit  on  AT^g* 

TmelfTin-  ln(l+t/r)  +  q"/h  +  (^)2Lq-/r^cpr  (2?) 

Tcorr’Tin  -  q"/h  +  (r+t)2Lq"AXJc  r  (28) 
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f-i  i  •  us  site  to  2  mm. 


14.  71-rt  fatigue-sputter  life  limits,  the  temperature 
limitS,  and  t>.e  minimum  tube  radius  limit  constitute  the 
thermo-mechanical  design  window  constraints. 

15.  For  flat  plates  with  circular  internal  channels, 
the  above  procedure  is  also  applicable,  except  that  the 
Quantity  (t+r)/r  does  not  ✓ary  as  the  plasma-facing  thick¬ 
ness  varies. 

16.  Protective  tile  which  can  expand  axially  but 
cannot  bow,  and  which  is  not  bonded  to  the  suostrate  is 
subject  to  the  same  tnermo-mechanical  constraints  as  above. 
For  the  tile  mechanically  bonded  to  the  substrate,  the 
analysis  of  equations  33  to  40  provides  the  expected  strain 
range  for  the  material.  Thickness  may  be  varied  until  the 
expected  sputtering  life  just  equals  the  fatigue  life  cor¬ 
responding  to  the  strain  range  at  that  thickness.  This 
provides  the  optimum  life  and  thickness  for  the  protective 
tile  bonded  to  the  substrate. 

1?.  Material  properties  are  input  for  average  temper¬ 
atures.  except  that  mechanical  limits  are  chosen  to  be  con¬ 
servative  over  the  temperature  range  of  interest.  Fluid 
properties  are  evaluated  at  the  bulk  or  wall  temperatures 
using  the  fits  of  empirical  data  in  Ref.  23. 


